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INTRODUCTION 


During  many  years  the  intensive  theoretical  and  experimental  investigations  of  a 
boundary  layers  existing  near  the  bodies  streamlined  by  the  flow  of  gas  or  fluid  were  studied. 
The  development  of  a  modern  aviation  and  especially  the  supersonic  one  requires  effective 
means  for  the  control  of  characteristics  of  flow  near  the  aeroplane’s  surface,  boundary  layer  heat 
and  mass  transfer  in  these  regions,  flow  separation  and  boundary  layer  friction  drag  reduction. 
With  the  purpose  of  improvement  of  the  aerodynamic  characteristics  of  flight  vehicles  the  set  of 
the  methods  of  an  external  influence  on  a  boundary  layer  was  offered.  However  existing  methods 
of  an  application  of  energy  to  a  boundary  layer  are  ineffective  or  economically  unprofitable. 
Therefore  the  development  of  new  ways  of  influence  on  the  characteristics  of  gas  flow  near  a 
flight  vehicle  exists. 

Recently  the  interest  to  use  for  various  applied  tasks  of  low-temperature  plasma  has 
revived  again.  So,  for  example,  a  new  direction  of  aerodynamics  -  so-called  plasma  supersonic 
aerodynamics  was  arisen  [1-4].  In  this  case  the  various  type  of  the  gas  discharges  are  used  with 
the  purpose  of  influence  on  the  characteristics  of  gas  flow  near  a  surface  of  the  flying  bodies. 
However  a  physics  of  the  discharge  in  supersonic  flow  of  gas  until  the  present  time  is  in  a  phase 
of  development.  There  are  many  unsolved  questions,  among  which  a  problem  of  a  gas 
breakdown  in  supersonic  flow,  a  creation  and  maintenance  of  the  stable  discharge  in  supersonic 
flow  of  air,  an  influence  of  flow  on  parameters  of  plasma  of  the  gas  discharge  and  an  influence 
of  the  discharge  on  the  characteristics  of  supersonic  flow  are  dominant  problems. 

The  first  laboratory  experiments  [1-3]  have  shown  an  opportunity  of  a  drag  reduction  at 
creation  of  the  discharges  of  direct  and  alternative  currents  before  a  body  streamlined  supersonic 
airflow.  However,  the  electrode  discharges  in  flow  are  unstable  and  spatially  non-uniform.  Such 
discharges  result  in  strong  erosion  of  electrodes  and  model  surface  and  reliably  are  not 
reproduced  in  various  realizations  [1]. 

A  task  of  search  of  optimum  ways  of  creation  of  non-equilibrium  plasma  in  supersonic 
flows  was  appeared.  One  of  such  ways,  offered  in  our  laboratory,  is  the  new  version  of  a  surface 
microwave  discharge,  namely,  the  microwave  discharge  on  an  external  surface  of  the  dielectric 
bodies,  streamlined  supersonic  flow  of  air.  In  this  case  the  economically  efficient  influence  on 
the  delay  of  laminar-turbulent  transition  and  on  the  reduction  of  a  turbulent  friction  coefficient  at 
the  external  action  on  a  turbulent  boundary  layer  can  take  place. 


The  Project  "Investigation  of  discharges  influence  on  boundary  layer  in  supersonic 
airflow"  is  devoted  to  the  development  of  methods  of  creation  of  plasma  of  different  types  of  the 
gas  discharges  near  the  surface  of  bodies  and  in  the  boundary  layers.  The  main  objective  is  to 
investigate  the  physics  of  non-equilibrium  plasma  of  the  surface  microwave  discharge  in 
supersonic  airflow. 


CHAPTER  I 


REVIEW  OF  THE  LITERATURE 

In  the  recent  years  the  task  of  improvement  of  the  aerodynamic  characteristics  of  flight 
vehicles  by  the  new  methods  is  an  actual  problem.  For  this  purpose  it  is  necessary  to  lower  a 
drag  reduction  and  to  reduce  a  friction  coefficient  in  a  boundary  layer. 

The  reduction  of  a  drag  can  be  achieved  by  creating  an  area  of  heated  gas  before  a  flight 
vehicle.  Such  area  can  be  formed  as  a  result  of  absorption  of  electromagnetic  energy  of  the 
different  frequency  ranges  or  by  plasma  creating  in  this  area.  For  this  purpose  in  some 
laboratories  it  is  offered  to  use  the  discharges  of  different  types. 

With  the  purpose  of  influence  on  the  characteristics  of  boundary-layer  flow  near  a  body, 
streamline  supersonic  gas  flow,  such  methods  as  heating  of  a  body  surface,  action  on  a  boundary 
layer  of  acoustic  oscillations,  application  of  vibrators,  suction  of  a  boundary  layer  etc.  are 
considered. 


1.1.  Frontal  drag  reduction 

The  system  of  shock  waves  existing  near  a  body  moving  in  the  Earth’s  atmosphere  with 
supersonic  velocity  causes  a  continuos  acceleration  of  air  masses  in  direction  of  the  flight.  It 
leads  to  the  frontal  drag  of  the  supersonic  aeroplane.  In  paper  [4]  it  is  offered  to  minimize  the 
acceleration  of  air  in  a  direction  of  a  motion  of  a  flight  vehicle.  For  it  the  airflow  should  be 
subsonic  everywhere.  For  elimination  of  frontal  shock  waves  it  is  enough  to  create  before  a  nose 
of  a  flight  vehicle  the  area  with  the  heightened  values  of  the  local  velocity  of  a  sound  vs,  for 
example,  to  heat  air  with  the  help  of  absorption  of  beam  energy.  It  is  important,  that  the  gas 
heating  should  not  be  bound  up  to  an  injection  of  a  preheated  continuous  medium  -  plasma  or 
gas:  otherwise  in  a  direction  of  "?"  on  boundary  with  cold  gas  the  local  velocity  of  a  sound  in 
cold  gas  will  be  exceeded,  the  shock  wave  will  be  formed  and  the  radical  drag  reduction  will  fail. 

It  is  necessary  to  generate  such  area  in  air  even  before  approach  of  a  flight  vehicle.  If  the 
shape  of  this  area  is  elongated  enough  on  "T  direction,  the  gas  expanding  at  a  heating  has  (in  a 
laboratory  frame)  practically  only  radial  component  of  velocity  vr.  This  velocity  depends  on  the 
radial  co-ordinate  and  the  time-space  allocation  of  an  energy  release  on  cross  section  at  co- 


ordinate  "?",  but  does  not  depend  on  velocity  of  a  flight  vehicle  v.  Apparently,  the  velocity  vr 
can  be  enough  small  in  comparison  not  only  with  v,  but  also  with  vs,  i.e.  it  is  possible  ensure  its 
subsonic  expansion  without  formation  of  strong  shock  waves.  If  the  area  of  an  energy  release 
insufficiently  elongates  (for  example,  quasi-spherical  zone)  it  is  streamlined  by  cold  gas  almost 
as  a  solid  body  and  the  shock  wave  is  formed. 

In  order  that  the  motion  of  a  flight  vehicle  in  an  atmosphere  happened  without  formation 
of  frontal  shock  waves  two  requirements  should  be  executed. 

First, 


where  k  is  Boltzmann  constant,  T  is  temperature,  m  is  mean  molecular  mass,  y  is  adiabatic  index. 

Second,  the  final  radial  dimension  of  the  gas-plasma  channel  rp  should  be  more  then  the 
radial  size  ra  of  a  nose  of  a  flight  vehicle. 

Let's  consider  the  problem  on  creation  of  the  heated  area.  It  can  be  generated  as  a  result 
of  absorption  of  electromagnetic  energy  of  the  different  frequency  bands  or  at  creation  of  plasma 
of  a  gas  discharge.  One  or  two  electrodes  are  required  for  formation  of  the  direct  current,  pulsed- 
periodic  or  high-frequency  discharges.  In  this  case  before  a  body  it  is  necessary  to  dispose 
electrodes  so  that  to  ensure  the  greatest  frontal  drag  reduction  at  a  relative  supersonic  motion  of  a 
body. 

In  paper  [5]  the  model  was  placed  along  a  wind  tunnel  axis,  the  cathode  was  placed  at  the 
front  of  model,  and  anode  was  behind  it.  In  experiment  two  types  of  models  (with  spherical  and 
conical  frontal  parts)  were  used.  The  anode  consisted  of  8  electrodes,  each  of  which  was 
included  to  the  power  source  through  the  ballast  resistance.  The  plasma  was  created  on  the 
frontal  surface  of  model.  Thus  the  frontal  drag  reduction  up  to  15  %  for  model  to  a  spherical 
surface  and  up  to  10  %  with  conical  one  took  place. 

The  diminution  of  a  frontal  drag  reduction  can  be  explained  by  two  mechanisms.  At  the 
considerable  gas  heating  the  translational  temperature  of  molecules  has  increased  up  to  1000- 
2000  ?.  Therefore  for  some  modes  the  drifted  plasma  has  a  subsonic  speed,  that  gives  in  an 
uniform  distribution  of  pressure  in  plasma,  and  to  absence  of  shock  waves  in  plasma  area.  The 
Reynolds’s  numbers  are  much  lower  than  for  the  same  stream  without  the  discharge.  The 
creation  of  the  discharge  probably  hinders  with  laminar-turbulent  transition  and,  thus,  reduces  to 
a  viscous  component. 

Attempt  to  investigate  a  degree  of  influence  of  one  of  these  mechanisms  on  a  frontal  drag 
reduction  was  conducted  and  enunciated  in  paper  [6].  During  experiment  the  efforts  were  guided 


on  separation  of  plasma  and  thermal  effects.  But  the  obtained  data  have  not  yielded  results  due  to 
the  large  errors. 

In  paper  [7]  the  chemical  reactions  were  taken  into  account  at  examination  of  plasma 
effects  in  non-equilibrium  plasma.  The  detailed  experimental  and  theoretical  examination  of 
plasmachemical  reactions  which  are  taken  place  at  a  microwave  discharge  in  air  enunciated  in 
the  monograph  [8].  In  these  papers  the  conditions  of  realization  of  experiment  are  close  to 
requirements  of  the  plasma  aerodynamic  experiments. 

At  usage  of  the  microwave  or  the  laser  radiations  they  spoke  about  an  electrodeless 
discharge.  The  adaptability  to  manufacture  of  energy  transportation  and  its  input  into  the 
discharge  chamber,  the  efficiency  of  absorption  of  electromagnetic  energy  in  the  discharge 
plasma,  the  absence  of  plasma  instability  initiated  by  electrodes,  the  long-lived  maintenance  of 
cleanness  of  plasma  because  of  absence  of  contact  with  electrodes  and  walls  of  the  discharge 
chamber,  the  opportunity  of  creation  of  the  extended  homogeneous  discharge,  the  achievement 
of  the  high  specific  supplied  powers,  the  selectivity  at  redistribution  of  energy  from  electrons  in 
the  internal  degrees  of  freedom  of  molecules  are  the  features  of  a  microwave. 

The  application  of  a  microwave  energy  to  plasma  can  be  executed  as  follows  [9].  The 
dielectric  tube  from  a  material  transparent  for  a  microwave  (for  example,  quartz)  is  interposed 
into  a  waveguide  of  the  rectangular  section  perpendicularly  to  its  broad  wall.  The  vector  of  an 
electric  field  in  a  wave  such  as  ?oi  is  parallel  to  the  narrow  wall  of  a  waveguide,  the  electric 
field  ?  in  this  direction  is  constant,  and  in  perpendicular  direction  it  changes  under  the  sine  law. 
The  discharge  is  formed  in  the  zone  of  cross  of  a  tube  with  a  waveguide.  The  plasma  is 
supported  at  the  expense  of  a  dissipation  of  energy  of  a  running  electromagnetic  wave. 

However  the  peculiar  properties  of  a  microwave  allows  to  conduct  experiments  in  free 
space  in  absence  of  a  contact  of  the  discharge  with  walls  of  the  chamber  and  electrodes.  In  [10] 
the  discharge  formed  in  the  resonator  in  which  the  standing  electromagnetic  waves  such  as  ?oi 
were  excited.  The  discharge  was  formed  on  an  axis  in  the  zone  of  a  maximum  electric  field.  The 
plasma  was  drawn  out  along  vector  of  an  electric  field  and  looked  like  a  cord  at  large  powers. 

The  successes  in  development  of  an  electronics  of  ultrahigh  frequencies  have  resulted  in 
an  opportunity  of  creation  of  the  new  form  of  the  discharge  -  the  electrodeless  microwave 
discharge  in  the  focused  beam  of  electromagnetic  energy  in  a  free  space.  In  [11]  the  microwave 
discharge  was  observed  at  argon  pressure  p=100  torr.  The  focused  beam  was  formed  with  the 
help  of  a  mirror.  The  radiation  of  the  microwave  generator  working  in  a  continuous  mode  power 
W=2  kW  on  a  wavelength  A.=3  cm  was  used.  The  pulsed  microwave  discharge  in  the  focused 
beam  of  radiation  was  explored  in  [12]  too. 


The  microwave  discharge  arising  in  the  focal  zone  of  electromagnetic  energy  beam 
represents  a  complex  nonlinear  phenomenon  including  a  non-stationary  breakdown  of  gas, 
spreading  of  ionization  fronts  interacting  with  electromagnetic  radiation,  maintenance  of 
nonequilibrium  plasma,  heating  and  excitation  of  neutral  gas.  Under  action  of  the  powerful 
focused  beam  of  electromagnetic  waves  of  a  very  high  frequency,  provided  that  the  electric  field 
intensity  exceeds  threshold  value,  a  gas  breakdown  takes  place  in  the  focal  zone  of  the  discharge 
chamber.  A  microwave  discharge  in  the  focused  beam  was  investigated  in  detail  in  a  cycle  of 
paper  [8,  13-23]. 


1.2.  Reduction  of  boundary  friction 

The  laminarization  of  airflow  is  one  of  the  few  possibilities  of  reducing  of  a  friction 
resistance  and,  consequently,  improving  the  aerodynamic  quality  and  efficiency  of  flight 
vehicles.  For  subsonic  aircraft  the  friction  resistance  makes  up  till  80-90  %,  and  for  supersonic 
till  50-60  %  from  total  resistance.  It  is  known,  that  the  friction  coefficient  at  laminar  flow  almost 
on  the  order  is  less  than  at  turbulent  one.  From  here,  the  importance  of  a  problem  of  a 
laminarization  of  airflow  is  obvious  for  flight  vehicles.  In  recent  years,  much  attention  has  been 
given  to  the  study  of  new  possibilities  for  laminarization  of  flow,  i.e.  delay  of  transition  to  the 
turbulent  flow  mode  in  the  boundary  layer  up  to  the  highest  possible  Reynolds’s  numbers.  One 
of  the  unconventional  methods  of  controlling  the  laminar-turbulent  transition  [24]  is  the  heating 
of  a  small  area  of  the  surface  of  a  body  being  flown  around. 

It  has  been  established  in  [24]  that  the  heating  of  a  part  of  a  body  surface  in  the  vicinity 
of  the  leading  edge  helps  to  considerably  increase  the  extent  of  the  laminar  portion  of  boundary 
layer.  In  paper  [24]  was  shown,  that  the  heating  of  a  small  surface  area  in  the  vicinity  of  the 
leading  edge  to  temperatures  two  to  four  times  higher  than  the  forward  flow  temperatures  leads 
to  delay  of  transition  even  if  the  body  surface  is  heat-insulated  downstream  of  a  heating  area.  An 
example  is  provided  of  the  surface  temperature  distribution,  when  the  energy  savings  due  to 
reduction  of  friction  exceed  the  energy  spent  to  heat  the  surface  by  a  factor  of  three. 

In  [25]  the  results  of  experimental  examination  of  an  opportunity  of  a  laminarization  of  a 
boundary  layer  by  the  method  of  suppression  of  Tollmien-Schlichting  waves,  existing  in  a 
boundary  layer  and  responsible  for  transition,  were  represented.  Thus,  the  artificial  created 
waves  of  perturbation  were  inlet  into  a  boundary  layer  in  an  antiphase  to  Tollmien-Schlichting 
waves.  The  opportunity  of  such  laminarization  was  shown.  The  experiments  were  conducted  on 
a  direct  wing  in  a  low  turbulent  subsonic  wind  tunnel.  The  boundary  layer,  perceiving  external 


fluctuations,  quenches  them  in  a  broad  band  of  frequencies  except  for  the  some  narrow  strip. 
Thus,  it  is  possible  to  receive  diminution  of  amplitude  of  the  most  dangerous  vibration  by 
introduction  in  a  boundary  layer  the  artificial  created  waves  of  perturbation.  In  such  a  way,  a 
retardation  of  development  of  a  dominant  wave  of  instability  takes  place  and,  as  result,  the  line 
of  transition  displaces  downstream. 

The  problem  of  transition  of  a  laminar  boundary  layer  in  a  turbulent  condition  during 
several  decades  causes  a  big  interest  of  the  researchers.  It  is  stipulated,  at  first,  by  necessity  of 
the  solution  of  practical  problems,  for  example,  the  problem  of  control  of  a  boundary  layer  with 
the  purpose  of  a  reduction  of  resistance  of  aircraft  and  calculation  of  its  aerodynamic 
characteristics.  Secondly,  the  study  of  process  of  appearance  of  a  turbulence  is  a  part  of  more 
general  fundamental  problem  of  the  description  of  a  turbulence. 

In  the  literature  the  conception  of  "the  transition"  is  understood,  as  a  rule,  concept  alone 
process  of  decay  of  a  laminar  mode  and  formation  of  a  turbulent  flow.  A  beginning  of  transition 
usually  connects  to  occurrence  of  rough  processes  in  the  form  of  turbulent  spots  and  low- 
frequency  pulsations  of  large  amplitude.  However  it  became  now  absolutely  understandable,  that 
the  long  sequence  of  physical  processes,  resulting  in  a  laminar  flows  destroying,  takes  the 
beginning  much  earlier,  namely,  in  transformation  of  external  perturbations  of  a  different  nature 
in  an  external  flow  in  waves  of  a  boundary  layer  or  in  their  generation  into  the  boundary  layer 
happening  on  unevenness’  of  a  streamline  surface  (ledges,  roughness  etc.). 

At  the  end  of  XIX  century  Reynolds  expressed  the  assumption  that  the  instability  of 
laminar  flow  is  the  reason  of  transition  of  laminar  flow  in  a  turbulent  condition.  In  other  words, 
the  waves  in  a  boundary  layer  amplify  and  it  is  the  reason  of  destruction  of  a  laminar  mode  of 
stream.  Much  later  Taylor  was  expressed  other  hypothesis,  according  to  which  the  transition  is 
caused  by  pulsation’s  of  an  external  flow  resulting  in  to  local  boundary  layer  separation  and  its 
turbulization. 

The  theory  of  magnetohydrodynamic  stability  started  by  Relay’s  transactionses  on  a 
boundary  of  ?  ?  century.  The  first  calculations  of  boundary  layer  stability  were  carried  out  by 
Tollmien  and  Schlichting  much  later  at  the  end  of  20  -  at  the  beginning  of  30  years.  However 
right  up  to  40  years  Tailor’s  hypothesis  was  preferred  in  a  stability  theory.  Only  Schubauer’s  and 
Skramstad’s  experiments  carried  out  at  the  beginning  of  40  years  have  confirmed  validity  of  the 
concept  of  instability.  Schubauer  and  Skramstad  have  found  out  both  oscillation  of  a  boundary 
layer  and  their  determining  role  during  destruction  of  a  laminar  mode.  Now  it  is  absolutely  clear 
that  in  case  of  small  intensity  of  all  external  perturbations  and  smooth  surface  of  a  plate  the 


transition  to  a  turbulence  in  a  boundary  layer  happens  owing  to  development  of  instability  of  an 
initial  laminar  flow. 

The  process  of  transition  of  a  laminar  boundary  layer  in  a  turbulent  condition  at  small 
intensity  of  external  perturbations  consists  of  three  conditionally  parted  stages:  the  generation  of 
waves  of  a  boundary  layer,  their  amplification  under  the  laws  of  the  linear  theory  and  non-linear 
destruction  of  a  laminar  mode  of  flow.  Each  stage  of  this  sequence  corresponds  the  typical  area 
in  space  in  accordance  with  increase  of  distance  from  a  leading  edge  of  model.  The  described 
sequence  of  stages  of  transition  is  diagrammatically  shown  in  a  Fig.  1.1  [26]. 

In  [27]  the  opportunities  of  improvement  of  the  aerodynamic  characteristics  of  a  plate  are 
investigated  at  the  expense  of  a  combined  laminarization  including  the  optimization  of  pressure 
distribution  and  temperature  of  a  surface.  The  class  of  the  temperature  distribution  of  a  surface 
corresponding  to  a  heating  of  a  leading  edge  and  cooling  of  a  remaining  part  of  a  surface  up  to 
the  recovery  temperature  of  an  incident  flow  was  surveyed.  It  was  shown,  that  thus  it  was 
possible:  i)  to  increase  a  lift  coefficient  of  a  plate  at  invariable  resistance  and  ii)  to  lower 
resistance  and  to  increase  a  critical  Mach  number  at  invariable  lift.  In  [27]  was  shown,  that  the 
energy  expedient  reduction  of  resistance  can  achieve  by  a  heating  of  a  part  of  a  wing  surface  in  a 
combination  to  other  methods  of  flow  laminarization. 

The  behaviour  of  a  separated  flow  at  a  motion  of  flight  vehicles  is  considered  in  [28].  The 
aerodynamic  characteristics  of  flight  vehicles  are  determined  by  features  of  their  streamlining, 
that  is  relative  motion  of  a  surrounding  medium.  The  flow  separation  is  one  of  the  fundamental 
phenomena  of  a  mechanics  of  a  fluid  and  gas.  It  consists  that  under  certain  conditions  (for 
example,  unfavourable  gradient  of  pressure)  the  medium  in  vicinity  of  a  streamlined  body  ceases 
to  move  along  its  surface  and  starts  to  moves  away  from  it.  In  result  the  area  of  the  separation 
flow  will  be  formed.  The  flow  separation  renders  influence  on  a  manoeuvrability,  controllability 
and  aerodynamic  efficiency  of  a  vehicle.  The  behaviour  of  the  separation  flow  differs  for 
different  aeroplane  speeds,  wing  shapes  and  attack  angles.  At  a  separation  of  a  stationary  laminar 
stream  the  separated  flow  can  again  join  a  wing  surface.  Thus,  the  local  separation  zone  arises.  It 
has  the  small  sizes  as  contrasted  to  the  wing  sizes  (Fig.  1.2,  a).  In  this  case  separated  area  on  a 
wing  exists  as  a  narrow  strip  elongate  along  its  total  wing-span  The  position  of  a  separated  zone 
on  a  wing  depends  on  an  angle  of  attack.  At  its  increasing  the  separation  zone  is  shifted  closer  to 
a  leading  edge  of  a  wing.  One  more  separation  of  a  non- stationary  turbulent  flow  is  appeared  on 
back  edge  of  a  wing  (Fig.  1.2,  b).  At  a  turbulent  separation  there  is  no  repeated  connection  of  a 
separated  flow  to  a  wing  surface  and  zone  of  a  separation  includes  all  area  of  flow  from  a  line  of 
a  separation  up  to  a  back  edge  of  a  wing.  A  flow  separation  from  a  leading  edge  arises  at  further 
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Fig.  1.1.  The  diagram  of  main  stages  of  transition  process  in  a  boundary  layer  [26] 


Fig.  1.2.  The  diagrams  of  the  flow  separation  at  different  values  of  an  attack  angle  of  a  wing  [28]. 


Fig.  1.3.  The  scheme  of  flow  in  vicinity  of  a  wing  surface  at  the  excitation  of  the  spatially 
homogeneous  acoustic  field  and  at  the  focused  acoustic  field  [28]. 


magnification  of  an  angle  of  attack  and  reaching  it  of  critical  value.  It  is  called  global  or 
complete  separation  (Fig.  1.2,  c).  The  complete  separation  is  accompanied  by  sharp  decrease  of 
lift  of  a  wing,  by  magnification  of  its  drag  and  gives  in  unfavourable  consequences. 

For  the  first  time  the  preventing  of  a  complete  separation  from  a  leading  edge  of  a  wing 
at  excitation  of  acoustic  oscillations  was  fulfilled  in  [29]  in  representation  about  two-dimension 
flow.  The  scheme  of  flow  in  vicinity  of  a  wing  surface  are  represented  in  Fig.  1.3.  The  excitation 
of  the  spatially  homogeneous  acoustic  field  from  a  remote  source  of  oscillations  causes  the 
connection  of  the  separated  flow.  In  result  the  complete  separation  from  a  leading  edge  is 
replaced  by  a  separation  of  a  turbulent  flow  downstream  (Fig.  1.3,  a,b).  Thus,  the  streamlining  of 
the  wing  remains  to  symmetric  concerning  its  central  section.  If  the  wing  is  in  the  focused 
acoustic  field  the  behaviour  of  the  vortexes  in  a  zone  of  a  separation  depends  on  what  area  of 
flow  is  put  to  the  test.  If  the  oscillations  focused  on  the  edge  of  model  they  have  an  effect  in  a 
zone  of  one  of  vortexes  and  the  influences  do  not  render  in  the  area  of  other  vortex  (Fig.  1.3,  c,d). 
From  here  it  is  possible  to  make  a  conclusion  that  the  spatial  structure  of  flow  at  acoustical 
action  and  the  result  of  control  of  a  complete  separation  depend  on  transport  of  gas  in  a 
transverse  direction,  in  spite  of  the  fact  that  the  acoustical  excitation  of  oscillations  of  the 
separated  layer  realizes  near  to  a  leading  edge  of  a  wing  and  the  large-scale  turbulent  motion  is 
formed  downstream. 

In  paper  [30]  the  experimental  study  of  a  susceptibility  of  a  boundary  layer  is 
investigated  at  modelling  of  process  of  excitation  of  own  oscillation  of  a  boundary  layer  by 
external,  artificial  caused,  acoustical  perturbations.  The  discharge-boundary  layer  system  was 
used  as  a  source  of  the  determined  perturbations.  The  surface  electric  discharge  in  a  laminar 
boundary  layer  of  first  plate  causes  the  perturbations  which  amplifies  in  some  times  at  a  motion 
downstream.  This  process  is  accompanied  by  acoustical  radiation  in  an  external  stream.  The 
distribution  of  the  perturbation  amplitude,  radiated  in  a  free  stream,  depends  on  longitudinal  co¬ 
ordinate  "?"  and  has  two  maximums.  The  first,  rather  narrow  maximum,  is  caused  by  electric 
discharge  and  is  spread  on  Mach  line  from  it.  The  second  maximum  corresponds  to  radiation 
from  a  transitional  zone  and  is  spread  on  a  line  by  close  to  Mach  line.  The  frequency  of  radiation 
coincides  frequency  of  an  electrical  arc  ignition.  When  the  second  plate  moves  in  a  field  of 
perturbations,  the  maximums  of  radiation  fall  on  different  parts  of  a  plate  and  cause  different 
response  of  a  boundary  layer  on  external  perturbations.  From  experiment  one  can  see  that  the 
most  intensive  transformation  of  external  acoustical  perturbations  to  own  oscillations  of  a 
boundary  layer  of  a  flat  plate  happens  in  neighbourhoods  of  a  leading  edge  of  a  plate,  sound 
branch  of  a  neutral  curve,  and  lower  branch  of  a  curve  of  neutral  stability. 


The  method  of  a  reduction  of  friction  based  on  a  heating  of  a  part  of  a  surface  of  a  flat 
plate  near  to  its  leading  edge  was  offered  in  [31].  In  this  case  the  diminution  of  a  total  resistance 
of  friction  on  a  plate  was  reached  due  to  increase  of  stability  of  gas  flow  heated  near  to  a  leading 
edge  and  further  travelling  above  more  cold  surface,  that  result  in  to  essential  magnification  of 
length  of  a  laminar  part  of  boundary  layer  flow.  However  at  large  enough  values  of  a  typical 
Reynolds’s  number  the  extent  of  a  turbulent  part  of  a  boundary  layer  considerably  exceeds  the 
length  of  a  laminar  part,  owing  to  what  in  this  case  further  diminution  of  a  friction  drag  is 
possible  only  by  a  reduction  of  the  turbulent  friction.  It  is  known  that  on  an  isothermal  surface, 
which  temperature  is  higher  than  the  flow  recovery  temperature,  the  friction  is  less  than  that  on 
an  adiabatic  surface.  However  uniform  heating  of  all  streamline  surface  is  bound  up  to 
considerable  technical  difficulties  and  in  this  case  the  energy  cost  of  a  boundary  layer  heating 
considerably  exceeds  the  profit  received  owing  to  a  friction  reduction. 

In  paper  [32]  the  results  of  numerical  investigations  of  influence  of  a  heating  of  a  local 
part  of  a  surface  of  a  flat  plate  on  a  flow  in  a  developed  turbulent  boundary  layer  are  submitted. 
The  surface  was  assumed  to  be  heat-insulated  downstream  of  the  heated  portion.  It  is  shown  that 
at  the  beginning  of  the  heating  portion  the  friction  is  increased  and  then  the  friction  coefficient 
decreases,  and  its  minimum  values  are  reached  on  the  heat-insulated  part  of  a  surface  following 
behind  the  heated  portion.  The  reduction  of  coefficient  of  turbulent  friction  as  contrasted  to  its 
value  on  an  isothermal  adiabatic  surface  is  spread  to  considerable  distance  downstream  from  a 
heated  portion.  Due  to  "long  memory  effect"  of  a  turbulent  boundary  layer  the  considerable 
friction  reduction  is  reached  at  local  heat  input  in  comparison  with  the  uniform  heating  of  all 
streamline  surface.  At  presence  of  several  following  one  after  another  the  heating  portion  this 
"long  memory  effect"  intensifies  and  more  essential  diminution  of  value  of  the  turbulent  friction 
can  be  obtained.  At  identical  total  energy  consumption  the  reduction  of  an  integrated  friction 
coefficient  almost  twice  is  more  than  in  case  of  an  uniformly  distributed  on  all  surface  energy 
supplied  to  gas  in  a  boundary  layer. 

The  article  [33]  is  development  of  the  paper  [32].  In  [33]  the  calculations  are  performed 
of  a  supersonic  boundary  layer  on  a  flat  plate  with  one  and  five  zones  of  local  heating  of  its 
surface.  Within  the  concept  of  "displacement  body",  allowance  is  made  for  the  interaction 
between  the  boundary  layer  and  external  nonviscous  flow.  It  is  shown  that  the  inclusion  of 
viscous- nonviscous  interaction  leads  to  significant  variation  in  the  distribution  of  the  local 
friction  coefficient. 


It  is  known  that  the  friction  resistance  decreases  under  conditions  of  heating  of  the 
surface  being  flown  about  [34,  35].  Based  on  the  results  of  calculation  of  a  compressible 
turbulent  boundary  layer  using  the  algebraic  model  of  turbulence,  it  was  shown  in  [36]  that  the 
efficiency  of  the  thermal  method  of  reducing  turbulent  friction  may  be  raised  considerably  owing 
to  localization  of  the  region  of  heat  input  to  the  boundary  layer.  In  [37]  the  theoretical  studies 
were  made  into  the  effect  of  the  thermal  properties  of  the  material  of  the  surface  being  flown 
about  on  the  reduction  of  turbulent  friction  under  conditions  of  local  heating  of  a  flat  plate. 

The  localization  of  regions  of  heat  input  to  gas  under  conditions  of  fairly  intense  heating 
leads  to  the  emergence  of  appreciable  gradients  of  the  displacement  thickness  of  boundary  layer. 
In  this  case,  the  interaction  between  external  nonviscous  flow  and  the  boundary  layer  results, 
even  under  conditions  of  flow  past  a  flat  surface,  in  the  emergence  of  a  longitudinal  pressure 
gradient.  The  presence  of  a  pressure  gradient  obviously  leads  to  a  variation  of  the  characteristics 
of  the  turbulent  boundary  layer.  In  view  of  this,  it  becomes  necessary  to  include  the  viscous- 
nonviscous  interaction  and  assess  its  effect  on  the  parameters  of  the  boundary  layer  in  the 
presence  of  localized  regions  of  heat  input. 

The  authors  of  the  article  [33]  treat  the  flow  past  a  plate  with  the  Mach  number  in  the 
incident  flow  =  2  and  the  velocity,  density,  temperature,  and  dynamic  viscosity  u^  ,pl  ,T* 
and  p^,  respectively.  The  initial  section  is  identified  in  the  region  of  the  developed  boundary 
layer,  which  is  characterized  by  the  value  of  Reynolds  number  Re0  =  u^,  0O  /pi  =4000  , 

calculated  by  the  parameters  of  incident  flow  and  momentum  thickness  in  this  section  0* 

Located  at  the  distance  L  downstream  of  the  above-identified  section  is  the  beginning  of 
the  heated  portion  of  the  plate  surface,  which  has  the  length  h  .  Two  options  will  be  treated 
below,  namely,  (a)  a  single  heated  portion,  and  (b)  five  series-arranged  heated  surface  portions  of 
the  same  length  h*  with  the  constant  distance  H  between  the  beginnings  of  two  neighboring 
portions.  Outside  of  the  heated  portions,  in  which  the  heat  flux  distribution  is  preassigned,  the 
plate  surface  is  assumed  to  be  adiabatic. 

The  simulation  of  the  compressible  turbulent  boundary  layer  is  performed  using  the  two- 
parameter  k-£  model  of  turbulence  [38].  For  numerical  solution  of  the  system  of  equations  of  the 
turbulent  boundary  layer  in  view  of  viscous- nonviscous  interaction,  in  [33]  dimensionless 
variables  (without  the  superscript  *)  are  introduced,  defined  by  the  following  relations: 
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Here,  S(x)  is  the  dimensionless  displacement  thickness  of  the  boundary  layer,  ue(x)  is 


the  velocity  on  its  external  boundary,  and  k*  and  £  ,  respectively,  denote  the  kinetic  energy  of 
turbulence  and  its  dissipation  rate.  Relations  (1.2)  further  include  the  characteristic  linear  size  1 
and  characteristic  Reynolds  number  Re,  which  may  be  preassigned  arbitrarily,  because  the 
parameters  of  the  developed  boundary  layer  in  the  initial  section  are  defined  by  the  value  of  the 
Reynolds  number. 

The  values  of  1  and  Re,  preassigned  in  the  calculations,  define  the  value  of  the 
momentum  thickness  in  the  initial  section  according  to  the  expression  9*  =  1  Re0/Re .  In  the 


/r 

calculations,  the  values  of  1  =  lm  ,  Re=2-10  and,  consequently,  90  =  9.992m  were  preassigned. 

In  order  to  include  the  interaction  between  the  turbulent  boundary  layer  and  external 
supersonic  flow,  in  [33]  the  concept  of  "displacement  body"  is  employed,  according  to  which  the 
effect  of  the  boundary  layer  on  external  flow  is  equivalent  to  nonviscous  flow  past  the  effective 
surface  formed  by  adding  the  displacement  thickness  of  boundary  layer  to  the  real  surface  of  the 
body. 

In  the  calculations,  constant  values  were  preassigned  for  the  distance  from  the  initial 
section  x0  to  the  beginning  of  the  first  heated  portion  L  =  250 9(l  =0.5  m  and  for  the  length  of 
each  portion  h  =  20  9 0  =  0.04  m.  The  distribution  of  heat  flux  qw  in  the  heated  portions  was 
preassigned  in  the  sine  form.  In  the  presence  of  one  heated  portion,  the  total  amount  of  heat, 
equal  to  the  integral  of  qw  over  the  portion  length,  was  Qw  =  1 0  4  (option  1)  and  2  x  1 0  4 
(option  2).  When  five  heated  portions  were  treated,  the  distance  between  the  beginnings  of  two 
adjacent  portions  was  H  =  199  90  =9.2  m,  and  the  amount  of  heat  supplied  to  every  one  of  them 
was  Qw  =19  4 . 


Fig.  1.4  gives  the  distribution  of  the  displacement  thickness  of  boundary  layer  and  of  the 
velocity  gradient  on  its  external  boundary  for  a  single  heated  portion  for  two  above-identified 
values  of  total  thermal  power  supplied  to  the  boundary  layer  [33].  The  dashed  curves  indicate, 


Fig. 1.4.  The  distribution  of  (1,  2)  the  displacement  thickness  of  boundary  layer  and  (3,  4)  of  the 
velocity  gradient  on  its  external  boundary:  (1,  3)  10‘4,  (2,  4)  2T0‘4. 

Dashed  curves  indicate  8  in  the  case  of  ignoring  the  viscous- nonviscous  interaction,  and 
continuous  curves  indicate  8  in  the  case  of  including  the  interaction  [33]. 


Fig.  1.5.  The  distribution  of  the  local  friction  coefficient  for  a  single  heated  portion: 

(1)  Qw^lO"4,  (2)  Qw=2T0‘4  [33]. 


for  comparison,  the  distribution  of  the  displacement  thickness  ignoring  the  viscous- nonviscous 
interaction.  Naturally,  the  displacement  thickness  in  the  heated  portion  and  behind  it  increases 
with  Qw.  The  interaction  between  the  boundary  layer  and  external  flow  results,  first,  in  a  smaller 

increase  of  the  displacement  thickness  than  in  the  case  of  ignoring  the  interaction  and,  second,  in 
the  upstream  propagation  of  the  effect  of  surface  heating  on  the  boundary  layer  thickness. 

As  a  result  of  viscous- nonviscous  interaction  before  the  boundary  layer,  a  fairly 
considerable  gradient  of  velocity  occurs  on  the  external  boundary  of  the  boundary  layer,  whose 
absolute  magnitude  is  the  greater,  the  higher  the  power  of  heat  input.  A  fairly  intensive 
acceleration  of  external  flow  is  observed  in  the  greater  part  of  the  heated  portion. 

The  longitudinal  gradient  of  pressure,  emerging  upon  inclusion  of  viscous-nonviscous 
interaction,  affects  considerably  the  local  friction,  as  is  shown  in  Fig.  1.5  [33]  which  gives  the 
distribution  of  the  local  friction  coefficient  for  two  values  of  Qw  defined  by  the  expression 


Peue  Re  8  ^  3y  Jw 

The  dot-and-dash  curve  in  Fig.  1.5  indicates  the  distribution  of  the  friction  coefficient  on  a 


fully  adiabatic  surface,  calculated  for  ue  =  1,  and  the  dashed  curves  indicate  the  distribution  of 


Cf  on  a  heated  surface  ignoring  the  interaction. 

The  deceleration  of  external  flow  before  the  heated  portion  is  accompanied  by  a 
considerable  decrease  of  the  friction  coefficient.  Due  to  acceleration  of  flow,  the  friction 
coefficient  throughout  almost  the  entire  heated  portion  is  higher  than  that  on  the  adiabatic 
surface.  This  increase  of  Cf  is  the  higher,  the  greater  the  heat  input.  The  inverse  dependence  is 


observed  when  the  interaction  is  ignored.  As  in  the  case  of  u  e  =  1 ,  the  maximum  reduction  of 
local  friction  is  attained  behind  the  heated  portion. 

The  interaction  between  the  boundary  layer  and  external  flow  has  a  marked  effect  both 
on  the  distribution  of  the  local  friction  coefficient  and  on  integral  friction.  The  total  reduction  of 
friction  due  to  local  surface  heating,  attained  at  x*,  =1.5m,  is  less  than  in  the  case  of  ignoring 
the  interaction  by  approximately  3%  with  Qw  =  10  4  and  by  10%  with  Qw  =2x10  4 . 


In  the  presence  of  several  heated  portions  with  one  and  the  same  value  of  Qw  the 
distribution  of  the  velocity  gradient  on  the  external  boundary  of  the  boundary  layer  in  the  region 
of  heat  input  is  almost  periodic,  as  is  seen  in  Fig.  1.6  [33].  A  gradual  downstream  decrease  is 
observed  of  the  maximum  (in  absolute  magnitude)  values  of  u'x.  The  dashed  curve  in  Fig.  1.6 
indicates  the  distribution  of  displacement  thickness,  calculated  ignoring  the  interaction. 


Fig. 1.6.  The  distribution  of  (1)  the  displacement  thickness  of  boundary  layer  and  (2)  velocity 
gradient  on  its  external  boundary  for  five  heated  portion  [33]. 


Fig.  1.7.  The  distribution  of  the  velocity  of  external  nonviscous  flow  [33]. 


Fig.  1.7  reflects  the  effect  of  the  boundary  layer  on  external  supersonic  flow  [33].  As  a 
result  of  periodic  deceleration-acceleration,  the  maximum  value  of  velocity  of  nonviscous  flow 
increases  slowly  downstream. 

Fig.  1.8  gives  the  distribution  of  the  local  friction  coefficient  both  ignoring  (dashed  curve) 
and  including  (continuous  curve)  the  interaction  for  five  heated  portions  [33].  The  dot-and-dash 
curve  indicates  the  distribution  of  the  friction  coefficient  on  the  adiabatic  surface.  No  qualitative 
changes  are  observed  in  the  behaviour  of  Cf  as  compared  with  the  case  of  a  single  heated 
portion. 

From  the  obtained  in  [33]  results  one  can  make  conclusions  that  the  interaction  between 
turbulent  boundary  layer  and  supersonic  flow  has  a  considerable  effect  on  both  local  and  integral 
friction.  This  effect  increases  with  the  power  of  heat  input  to  the  gas.  In  optimizing  the  thermal 
method  of  reducing  turbulent  friction,  based  on  local  heat  input  to  the  boundary  layer,  the 
viscous- nonviscous  interaction  must  be  included. 

1.3.  Surface  microwave  discharge 

One  of  the  ways  for  obtaining  of  a  stable  plasma  column  of  considerable  length  in  a  large 
range  of  gas  pressures  and  in  a  big  diapason  of  frequencies  is  use  of  the  surface  electromagnetic 
wave  of  the  microwave  region.  In  cycle  of  papers  [40-46]  the  problems  bound  up  with  creation 
of  surface  high-frequency  and  microwave  discharges  and  study  of  their  properties  are 
considered.  The  plasma  was  formed  in  a  discharge  tube,  the  energy  was  supplied  with  the  help  of 
a  wave  guide  through  the  special  device  -  surfatron. 

In  a  microwave  surfatron  the  discharge  represents  a  plasma  waveguide  on  which  the 
surface  wave  supporting  ionization  of  gas  is  distributed.  Thus  for  a  wide  range  of  frequencies  it 
is  possible  to  receive  a  plasma  column  of  significant  extent  and  forms  determined  by  a  dielectric 
wall  limiting  the  discharge.  A  high  electron  concentration  (exceeding  critical  density)  and  its 
appreciable  unhomogeneity  in  a  direction  of  a  surface  wave  distribution  are  main  feature  of  such 
discharge.  It  is  known  [41,42]  that  at  creation  of  the  high-frequency  and  microwave  discharges 
inside  dielectric  tube  filled  by  gas  at  low  pressure,  the  electromagnetic  energy  delivered  to 
system  is  transformed  to  a  surface  wave.  Thus,  there  is  a  self-sustaining  system  when  a  plasma 
medium  created  by  the  surface  wave  is  necessary  for  a  surface  wave  existence,  i.e.  the  presence 
of  plasma  is  a  necessary  condition  for  spreading  of  a  surface  wave.  The  surface  wave  is  travelled 
in  space  so  long  as  its  energy  is  sufficient  for  creation  of  plasma  with  an  electron  density  no  less 
then  a  critical  value.  Sufficient  microwave  absorption  for  plasma  generation  can  only  be 
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achieved  for  overcritical  electron  densities  ne>nc.  The  critical  density  is  determined  by  the 
angular  frequency  of  the  electromagnetic  field  GO  and  the  permittivity  of  the  free  space  £0  : 

nc  =  me£oar  /  e2  (1.4) 

(e  and  me  are  the  electron  charge  and  mass,  respectively).  This  way  of  plasma  production  and 
device  for  its  creation  refer  to  as  a  surfatron. 

The  feature  of  such  kind  of  the  surface  discharge  is  the  longitudinal  density  gradient 
dne/dz  of  electrons.  It  is  explained  by  a  consecutive  supply  of  a  microwave  energy  to  different 
segments  on  its  length.  At  not  too  small  gas  pressure  a  high  electron  concentration  and  its 
longitudinal  gradient  should  give  in  considerable  non-uniform  heating  of  gas.  This  heating  leads 
to  change  of  the  reduced  gas  density  and,  therefore,  of  kinetic  coefficients  in  plasma.  The  energy 
balance  in  the  discharge  changes  also. 

The  paper  [47]  is  specially  focused  on  large  area  microwave  plasma.  It  briefly  reviews 
recent  results  in  the  field.  Only  their  basic  principle  are  analyzed.  The  discharges  are  classified 
according  to  operating  principle:  large  diameter  plasma,  radiating  structures  and  based  on 
evanescent  waves  structures.  Basic  properties  and  operating  conditions  are  described.  Most  of 
these  structures  have  been  designed  as  plasma  reactors  dedicated  to  plasma  material  processing. 

Interest  for  plasma  processes  strongly  increased  in  the  90's  particularly  considering  the 
microwave  plasma.  A  very  large  variety  of  applications  is  investigated  as  materials  processing, 
plasma  processing,  plasma  chemistry,  depollution...  New  plasma  sources  have  been  designed  for 
use  in  such  applications.  At  low  (or  medium)  pressure,  microwave  plasma  produce  high  densities 
of  charged  particles  and  further  they  exhibit  exciting  qualities  for  the  manufacturer  as  versatility 
(without  contamination),  improved  throughput  and  reduced  processing  cost.  Surfaguides  have 
been  extensively  reported  [48-57]  but  only  a  few  recent  papers  dealt  with  the  production  of  large 
plasma  [58-60]. 

The  paper  [61]  discusses  basic  concepts  of  slotted  waveguide  microwave  plasma 
excitation  which  is  technically  very  reasonable  for  the  generation  of  long  extended  and  large  area 
plasmas.  Results  of  different  experiments  aimed  at  the  details  of  the  spatial  plasma  structure  are 
reported  for  an  improved  version  of  slotted  waveguide  plasma  excitation  which  allows  maximum 
suppression  of  plasma  homogeneity  limiting  effects. 

Microwave  plasma  generation  is  determined  by  the  specific  propagation  properties  of 
microwaves  and  by  their  strong  absorption  in  the  plasma.  This  is  of  special  importance  for 
plasmas  with  dimensions  comparing  to  or  exceeding  the  microwave  free  space  wavelength.  In  a 
comprehensive  comparison  of  radio  frequency  and  microwave  plasma  reactors  Moisan  et  al. 


pointed  out  that  regarding  a  constant  power  density  operation  it  is  always  better  to  operate  at 
high  frequencies  for  processes  with  a  threshold  energy  below  about  half  that  for  ionization  [62]. 

To  generate  a  plane  wavefront,  a  field  antenna  with  characteristic  dimensions  comparing 
to  the  microwave  window  dimension  has  to  be  operated  under  far  field  conditions.  Far  field 
conditions  are  characterized  by  the  relation 

r0>8(rs2A0)  (1.5) 

with  r0  is  distance  between  source  and  microwave  window,  rs  is  characteristic  source 
dimension,  and  is  the  wavelength  in  free  space  [63].  Following  this,  a  setup  is  necessary 
which  is  extremely  space  consuming.  Further,  shielding  against  leaky  waves  which  is  needed  for 
security  reasons  will  be  a  problem.  It  needs  additional  space  and  can  cause  unwanted 
waveguiding  effects  disturbing  plane  wave  propagation.  A  much  more  appropriate  solution  to 
spread  microwave  power  over  a  large  plasma  area  is  to  use  waveguiding  structures  which  are 
directly  coupled  to  the  plasma. 

A  very  convenient  method  often  used  method  for  power  spreading  are  slotted 
waveguides.  Slotted  waveguides  are  in  widespread  use  for  waveguide  couplers  and  antennas. 
The  basic  problem  of  spreading  microwave  power  over  large  plasma  areas  is  the  compensation 
of  power  losses  due  to  absorption.  For  the  linear  slow  wave  structures  a  small  tilt  between 
structure  and  vessel  was  used  to  improve  plasma  homogeneity.  For  slotted  waveguides  the  use  of 
resonance  effects  [64,  65]  and  tuning  of  the  power  coupling  efficiency  of  the  slots  [66,  67]  is 
reported. 

However,  in  any  case  slotted  waveguides  are  inhomogeneous  waveguides  by  principle. 
Characteristic  dimensions  of  the  slotted  waveguides  compare  to  characteristic  dimensions  of 
power  transfer  to  the  plasma.  The  question  is  to  which  extend  this  discreteness  of  the  coupling 
slots  will  limit  the  generation  of  homogeneous  plasmas. 

The  interesting  effect  of  a  slot  in  a  waveguide  wall  is  that  it  may  interrupt  surface 
currents  on  the  inner  side  of  the  waveguide  wall.  This  way  an  electric  field  can  be  formed  in  the 
slot  which  can  leak  out  or  radiate  to  the  surrounding.  The  level  of  the  radiated  power,  its 
radiation  pattern  and  its  feedback  on  the  waveguide  depend  on  size  and  location  of  the  slot  and 
on  the  type  of  transmission  mode  in  the  waveguide.  Here  only  a  brief  summary  necessary  for 
basic  understanding  is  given.  It  refers  to  the  case  of  slots  in  the  broad  plane  of  a  rectangular 
waveguide  with  Hu  -  transmission  mode  [68].  The  situation  for  a  general  inclined-displaced  slot 
is  illustrated  in  Fig.  1.9  [61]. 

A  single  narrow  rectangular  slot  acts  like  a  dipole  radiator.  When  the  slot  has  been  tuned 


centre  line 


Fig.  1.9.  Scheme  of  a  section  of  a  slotted  rectangular  waveguide 
showing  a  general  inclined-displaced  slot. 

On  the  inner  walls  the  wall  current 
distribution  for  H10  -  mode  traveling  waves  is  indicated  [61]. 


by  cutting  to  a  length  ls ,  close  to  one  half  of  the  free  space  wavelength  it  is  resonant  and 
presents  a  pure  resistance  to  the  wave.  Therefore,  it  has  no  influence  on  the  phase  velocity  of 
wave  propagation.  A  longitudinal  slot  (cp  =  0)  presents  a  shunt  load.  It  is  excited  by  the 
transverse  component  of  the  surface  current  in  the  inner  side  of  the  waveguide  wall.  This  means 
that  it  remains  unexcited  if  the  displacement  dd  to  the  center  line  of  the  broad  face  is  zero.  The 
transverse  component  jst  of  the  surface  current  increases  with  the  distance  dd  to  this  line 
(jst  -  cos  Tlx /a)  and,  therefore,  the  power  coupling  increases  (~  cos2  Tlx /a).  A  transverse  slot 
symmetrically  cut  in  the  broad  face  (cp  =  90°,dd  =  o)  is  excited  by  the  longitudinal  component  of 

the  surface  current.  It  presents  a  series  load  to  the  wave  in  the  waveguide.  It  may  be  rotated 
about  its  center  and  still  remains  a  series  load  while  the  coupled  power  varies  from  zero  (at 

cp  =  0 )  to  maximum  (at  cp  =  90° ). 

Since  resonant  slots  present  pure  resistances  they  have  no  influence  on  the  phase  velocity 
of  wave  propagation  in  the  waveguide.  This  is  the  reason  why  they  are  often  used  to  construct 
antenna  arrays.  If  the  single  slots  are  in  phase,  the  radiation  pattern  formed  by  the  superposition 
of  the  single  slot  dipole  pattern  has  a  maximum  in  the  direction  normal  to  the  broad  plane  of  the 
waveguide  can  be  obtained.  This  is  the  case  for  slots  which  are  arranged  at  a  distance  of  about 
one  half  of  the  wavelength  in  the  waveguide  A,h  /  2 .  If  an  even  number  of  slots  are  arranged  at 
distances  of  a  quarter  wavelength  /4  the  array  radiates  in  the  forward  direction  of  the  wave, 
basically.  Similarly,  directional  waveguide  couplers  can  be  constructed  on  the  basis  of  resonant 
slots  which  transmit  power  in  the  forward  direction  only. 

In  the  simplest  case  the  dimensions  of  slotted  waveguides  used  for  microwave  discharge 
excitation  are  similar  to  linear-array  slotted  waveguide  antennas  for  radiating  waves  [69].  The 
only  difference  is  that  a  discharge  vessel  with  a  microwave  transparent  wall  is  placed  parallel  to 
the  waveguide  at  a  certain  distance  d  ,  which  is  short  compared  to  the  far  field  condition.  This 

way,  the  antenna  becomes  a  nearly  perfect  transmission  line  applicator. 

Slotted  waveguide  microwave  applicators  are  intended  to  generate  large  area  long 
extended  plasmas.  Therefore,  large  microwave  windows  with  considerable  thicknesses  are 
necessary.  Especially,  if  these  windows  are  inserted  in  a  conducting  vessel  wall  they  can  act  as 
an  additional  waveguide  which  is  now  inserted  between  the  waveguide  face  and  the  plasma  and 
which  can  be  resonant. 

Another  design  problem  for  long  linear  microwave  field  applicators  and,  therefore,  also 
for  slotted  waveguide  microwave  applicators  is  to  get  a  long  extended  plasma  with  maximum 


homogeneity.  Special  care  has  to  be  taken  to  ensure  homogeneous  power  distribution  since 
power  absorption  along  the  waveguide  is  not  negligible.  In  this  context,  it  is  necessary  to 
subdivide  between  two  configurations  of  the  waveguiding,  power  distributing  field  applicators, 
the  one  which  carries  standing  waves  and  the  other  which  carries  travelling  waves.  The 
travelling-wave  structure  is  terminated  by  a  matched  load.  The  standing-wave  structure  is 
terminated  by  fully  reflecting,  usually  tunable,  loads.  Both  configurations  possess  advantages 
and  disadvantages. 

There  are  different  methods  to  improve  slotted  waveguide  applicator  performance.  It  is 
well  known,  that  external  static  magnetic  field  can  be  used  to  control  the  homogeneity  of 
microwave  plasma.  The  use  of  magnetic  fields,  however,  is  limited  to  lower  pressures  <0. 1  mbar. 
At  higher  pressures  the  electron  cyclotron  resonance  effect  is  suppressed  by  collisional  damping. 
Therefore,  other  means  to  influence  plasma  homogeneity  are  necessary. 

Yasuietal.  [65]  describe  a  standing-wave  applicator  where  an  adjustable  T-shaped  ridge 
is  inserted  in  the  slotted  waveguide.  It  is  used  to  control  the  electromagnetic  field  pattern  in  the 
waveguide. 

Another  method  to  improve  the  performance  of  standing-wave  applicators  is  the  use  of 
closed  loop  waveguides.  This  has  the  advantage  that  standing  wave  modes  with  very  regular 
power  distribution  to  the  single  slots  can  be  obtained.  In  [64,70]  a  setup  using  circular 
waveguides  with  circumferential  lengths  of  the  middle  line  of  the  cavity  of  about  integer 
multiples  of  the  waveguide  wavelength  is  described. 

The  examples  reported  in  [64,  67,  69,  70]  indicate  that  an  appropriate  arrangement  of  the 
gap  between  slotted  waveguide  face  and  discharge  vessel  surface  is  of  great  importance  for  the 
obtained  plasma  distribution.  It  is,  therefore,  strongly  recommended  to  use  microwave 
applicators  which  exhibit  well  defined  microwave  propagation  properties  in  this  gap.  A 
traveling-wave  applicator  where  this  is  realized  was  first  described  in  [66] . 

The  basic  problem  of  large  area  homogenous  microwave  plasma  generation  is  to  avoid 
inhomogeneous  power  coupling.  This  concerns  the  scheme  of  power  distribution  over  the  whole 
plasma  area  and  the  localization  of  power  absorption  in  the  plasma.  An  ideal  microwave 
applicator  should  generate  a  microwave  field  distribution  which  is  homogeneous  over  the  entire 
microwave  window.  It  was  the  aim  of  the  development  of  the  so  called  "planar  microwave 
plasma  applicator"  to  adapt  the  technical  advantages  of  slotted  waveguide  power  distribution  to 
this  requirement.  The  design  of  this  applicator  is  based  on  a  principle  of  distributed  and  tunable 
coupling  of  microwave  power  in  conjunction  with  the  use  of  an  additional  smoothing  waveguide 
which  supports  homogeneous  plasma  excitation.  In  Fig.  1.10  and  Fir.  1.11  the  technical 
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Fig. 1.10.  Schematic  description  of  a  planar  microwave  plasma  source 
with  a  modified  slotted  waveguide  microwave  field  applicator  consisting 
of  two  rectangular  waveguide  in  T-shape  configuration  [61]. 


Fig.1.11.  Scheme  of  the  side  view  of  a  long  planar  plasma  source 
with  the  modified  slotted  waveguide  applicator  described  in  Fig. 1.10  [61]. 


realization  of  this  principle  using  two  rectangular  waveguides  in  T-shape  configuration  is 
described  schematically.  A  modified  traveling- wave  slotted  waveguide  is  used  to  distribute 
microwave  power  over  the  desired  length.  Instead  of  rectangular  slots  coupling  holes  at  distances 
of  A,h  /  4  are  used.  These  holes  contain  movable  metallic  hooks  which  improve  the  efficiency 
and  the  tuning  range  of  power  coupling  considerably.  Since  the  hooks  are  polarized  elements, 
they  can  be  used  to  determine  the  orientation  of  the  coupled  microwave  field. 

Microwaves  coupled  to  this  additional,  so  called  "interface  waveguide"  can  propagate  in 
both  directions  to  the  ends  of  the  applicator.  These  ends  are  terminated  by  matched  loads  which 
are  located  above  the  ends  of  the  microwave  window  to  avoid  resonance  effects  of  the  window 
slot.  Following  this,  only  travelling  waves  exist  in  the  interface  waveguide.  This  way,  a 
smoothing  effect  is  obtained.  At  each  position  of  the  waveguide,  field  components  originating 
from  different  coupling  points  superpose  to  a  complex  wave  package  without  any  interferences. 
Therefore,  for  positions  between  the  coupling  points  the  effect  of  power  losses  due  to  the  strong 
absorption  in  the  plasma  can  be  reduced.  Further,  power  losses  of  the  wave  in  the  distributing 
waveguide  can  be  compensated  by  tuning  of  the  coupling  hooks.  This  way,  a  good  plasma 
homogeneity  over  the  whole  length  of  the  plasma  can  be  obtained. 

Finite-area  plane  microwave  plasma  source  has  been  designed  by  H.Sugai  and  coworkers 
[71].  A  pair  of  slots  formed  on  a  waveguide  base  plane  couples  the  microwave  power  (2.45  GHz 
-  1  kW)  to  the  plasma  as  shown  in  Fig.  1.12.  These  two  slots  (width  =  10  mm  and  length  =  61 
mm  =  half  a  wavelength)  are  separated  by  79.6  mm,  a  half  of  guide  wavelength.  Angle  of  slots 
with  respect  to  the  waveguide  axis  is  27°  (matching  of  the  radiation  resistance  with  the 
characteristic  resistance  of  the  guide).  Slots  plane  is  separated  from  the  cylindrical  plasma 
chamber  (diameter  =  22  cm)  by  a  17  mm  thick  quartz  window.  Plasma  is  created  in  argon  in  the 
range  20  to  350  Pa  (gas  flow  is  715  seem).  Axial  field  distribution  is  found  quasiconstant  from 
the  quartz  window  surface  (z=0)  at  low  density  over  8  mm  and  sharply  decreasing  as  the  electron 
density  becomes  higher  than  the  critical  one  (7.410  cm"  ).  The  measured  electron  temperature 
is  Te~2  eV  for  these  high  density  plasma.  As  these  surface  waves  are  evanescent  in  the  axial 
direction  and  propagate  in  the  r-0  directions  on  a  large  plasma  surface,  the  surface  wave 
excitation  from  small  slots  antennas  would  provide  a  new  method  for  production  of  large  "flat 
panef  plasma. 

The  SLAN  structure  proposed  by  J.Engemann  [72]  is  based  on  the  principle  of  a  slotted 
radiating  guide  whose  the  slots  are  perpendicular  to  the  waveguide  axis.  But  SLAN  system 
consists  of  an  annular  waveguide,  called  a  ring  resonator,  whose  slots  antennas  are  parallel  to  the 


Fig. 1.12.  Finite  area  planar  microwave  plasma  [71]. 
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Fig. 1.13.  SLAN  structure.  Chamber  diameter  =  16  cm  [72]. 


ring  axis  and  positioned  at  regular  intervals.  Basic  structure  (Fig. 1.13)  is  such  the  middle  line  of 
the  annular  waveguide  has  a  circumferential  length  of  five  waveguide  wavelengths.  Five 
coupling  slots  are  on  the  inner  side  of  the  nng  resonator  and  are  half  of  a  free-space  wavelength 
long.  Diameter  of  the  ring  is  chosen  to  develop  a  standing  TEio  wave  inside.  Plasma  tube  is 
positioned  inside  and  coaxially  to  the  ring.  A  coupling  probe  transfers  the  magnetron  power  to 
the  ring  and  its  distance  from  the  nearest  slot  is  a  quarter  of  waveguide  wavelength.  Hence,  the 
standing  wave  nodes  of  the  electric  field  are  at  the  slots  positions  and,  in  vacuum,  each  slot 
antenna  constructively  interfere  forming  in  the  plasma  chamber  a  TE51  mode  as  in  a  cylindrical 
cavity. 

The  all  structures  producing  large  diameter  plasma  exhibit  several  common  features. 
They  are  multimodes  structures,  nevertheless  the  plasma  is  produced  along  a  dominant  mode 
which  is  determined  by  the  applicator  geometrical  configuration,  the  microwave  power  and  the 
pressure.  The  spatial  distribution  of  the  microwave  power  is  never  homogeneous,  hence  the 
homogenization  of  the  plasma  strongly  depends  on  the  diffusion  efficiency.  As  a  consequence, 
increasing  the  plasma  diameter  requires  lower  operating  pressures.  In  all  structures,  an 
appropriate  magnetic  field  can  be  added  either  for  obtaining  the  ECR  conditions  or  only  a  plasma 
confinement. 

General  principle  of  the  radiating  structures  is  based  on  gas  ionization  by  a  radiated 
wave.  This  wave  is  produced  either  by  slots  or  by  antennas,  the  excitator  being  designed  in  such 
a  way  as  every  slot  antenna  is  radiating  the  same  power.  Further,  the  dimensions  of  the  plasma 
chamber  are  much  larger  than  the  wavelength  of  the  radiated  wave.  Hence,  on  the  opposite  of  the 
large  diameter  plasma,  no  eigenmodes  can  occur  and  the  microwave  field  distribution  does  not 
depend  on  the  chamber  geometry.  We  briefly  present  four  of  these  structures. 

Fig.  1.14  shows  the  principle  of  the  planar  microwave  plasma  source.  It  consists  of  two 
waveguides  arranged  parallel  to  each  other  in  a  T-shape  cross-section  configuration.  Waveguide 
1  is  a  guide  on  one  end  and  is  connected  with  the  microwave  generator  on  the  other  end  with  a 
matched  load.  Waveguide  2  is  connected  at  both  ends  to  matched  loads.  Coupling  of  microwave 
power  from  guide  1  to  guide  2  is  obtained  by  means  of  radiating  antennas  shaped  into  hooks. 
These  hooks  act  as  discrete  adjustable  coupling  elements  distributed  along  the  guides  and  enable 
the  ignition  of  an  homogeneous  plasma  provided  the  elements  yield  the  same  power  (radiated 
waves  from  every  element  overlap  to  an  integral  field  distribution  of  nearly  constant  amplitude 
over  the  whole  plasma  length).  Obviously,  the  lower  side  of  guide  2  is  replaced  by  a  quartz 
window  separating  it  from  the  plasma  vessel,  the  plasma  itself  acting  as  a  conducting  wall.  The 
structure  operates  at  2.45  Ghz,  in  the  pressure  range  1  Pa  to  10  kPa  with  argon,  hydrogen  plasma 
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Fig.  1.14.  Planar  microwave  plasma  source  [47]. 


Fig. 1.15.  Resonant  field  applicators  [47] 

(a)  array  of  centred,  inclined  series- slots  in  the  broad  wall  of  a  waveguide 
(b)  array  of  longitudinal  shunt-slots  placed  off  the  centre  line  of  the  broad  wall  of  a  waveguide. 


with  power  ranging  from  0.1  to  1  kW.  Gases  are  injected  near  the  quartz  window  with 
maximum  flows  of  103  seem.  Homogeneity  better  than  5%  has  been  obtained  over  30x100  cm 
plasma.  It  is  expected  that  scaling-up  to  200  cm  could  be  obtained. 

Fig.  1.15  shows  the  principle  of  the  resonant  field  applicators.  A  waveguide  terminated  by 
a  short-circuit  plane  establishes  a  standing  wave  within  it.  It  delivers  power  to  the  plasma  by 
means  of  slots  at  positions  closely  correlated  to  the  standing  wave  pattern.  The  configuration  of 
the  slots  ensures  that  the  phase  of  the  field  radiated  by  each  slot  is  the  same  yielding  a  broadside 
radiation  pattern.  The  slots  are  located  at  minima  or  maxima  of  the  electric  field  and  their 
dimensions  ensure  that  slots  susceptance  or  reactance  vanishes  at  the  operating  frequency  (so- 
called  resonant  slots).  Two  such  configurations  are  shown  in  Fig.  1.15.  The  length  of  each  slot  is 
one  half  of  the  free-space  wavelength  X  (hence  the  resonance)  and  the  slots  are  spaced  one  half 
of  the  waveguide  wavelength  X„.  In  Fig.  1.15  a,  the  applicator  aperture  consists  of  an  array  of 
centred,  inclined  series-slots  in  the  broad  wall  of  the  waveguide,  the  short-circuit  being  located 
XJ2  past  the  last  slot  in  the  array.  In  Fig.  1.15  b,  the  aperture  consists  of  longitudinal  shunt-slots 
placed  off  the  centre  line  of  the  broad  wall  of  the  guide,  the  short-circuit  being  located  XJ4  past 
the  last  slot  in  the  array.  More  details  could  be  found  in  [73,  74]. 

Gigatron  and  Duoplasmaline  structures  are  described  in  more  details  in  [75].  Gigatron  is 
an  example  of  structure  based  both  on  radiating  structures  principle  and  propagation  of 
evanescent  surface  waves.  As  shown  in  Fig.  1.16,  it  consists  of  an  inner  copper  rod  surrounded  by 
a  quartz  tube  filled  with  atmospheric  pressure  air.  The  rod  radiates  the  microwave  power  coming 
from  a  coaxial  line  (supplied  by  a  2.45  Ghz  magnetron)  the  centre  conductor  of  which  is  the  rod. 
Either  one  or  two  magnetron  sources  can  be  used  giving  the  name  of  "Gigatron"  or 
"Duoplasmaline".  The  radiated  wave  yields  an  evanescent  surface  wave  which  propagates  along 
the  dielectric  interface.  The  strong  absorption  of  the  wave  yields  the  gas  ionization  in  the 
surrounding  chamber.  As  for  surface  waves  produced  plasma,  the  plasma  length  depends  on  the 
microwave  power.  Therefore,  with  the  Gigatron,  the  plasma  density  decreases  along  the  tube 
whereas  this  axial  decrease  is  compensated  by  the  use  of  two  magnetron  sources  with  the 
Duoplasmaline. 

Two  types  of  structures  based  on  evanescent  waves  are  described  [76-78].  Surface  wave 
along  a  dielectric  line  applicator  is  schematically  shown  in  Fig.  1.17.  A  rectangular  microwave 
waveguide  is  ended  by  a  fluorocarbon  polymer  sheet  covering  an  aluminum  plate  and  they 
consist  the  dielectric  line.  The  wave  coming  from  the  guide  generates  a  surface  wave  which 
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Fig. 1.16.  Duoplasma  line  structure  [75]. 


Fig.  1.17.  Structure  of  the  surface  wave  dielectric  line  [76]. 


propagates  along  the  dielectric  line.  This  wave  couples  the  microwave  power  to  the  gas  chamber 
via  the  evanescent  field  along  the  glass  window  separating  the  line  from  the  plasma  chamber. 
Profiles  of  the  electric  field  in  the  chamber  exhibit  a  standing  wave  pattern  imposed  by  the 
chamber  walls.  Although  this  pattern  is  complex,  it  has  been  found  a  good  agreement  with 
computations  assuming,  in  a  first  approximation  that  the  chamber  is  a  rectangular  guide  where 
the  microwave  mode  is  a  complex  mode  consisting  of  TM4m  and  TMgm  without  plasma.  Almost 
the  all  microwave  power  is  absorbed  near  the  microwave  window  producing  a  180x300  mm“ 
plasma,  the  pressure  ranging  in  the  10  to  100  Pa. 

The  coaxial-type  open-ended  dielectric  cavity  structure,  shown  in  Fig.  1.1 8  is  described  in 
more  details  in  [79].  The  microwave  power  of  a  TEio  mode  within  a  rectangular  guide  is  coupled 
to  a  TEM  coaxial  waveguide  mode  in  the  same  way  as  in  rectangular  coaxial  transitions.  Coaxial 
guide  is  then  widened  and  funnel-shaped  upside-down.  The  bottom  of  this  part  is  opened  on  a 
coaxial-type  cavity  whose  typical  dimensions  of  inner  and  outer  diameters  199  and  254  mm 
respectively.  At  the  open  end  of  the  coaxial-type  cavity  is  placed  a  Pyrex  window  of  20  mm 
thickness  and  175  mm  diameter  in  contact  with  the  inner  conductor.  Finally,  the  cavity  gap  mns 
into  the  circumferential  side  of  the  Pyrex  window.  Accounting  for  the  cylindrical  TEM  mode 
propagating  in  the  coaxial  guide,  the  Pyrex  window  is  the  dielectric  line  in  which  microwave 
propagate  from  the  edge  towards  the  centre.  Therefore,  the  surface  wave  radiates  into  the 
chamber.  The  microwave  power  decreases  from  the  edge  towards  the  centre  but  the  waves 
interfere  with  another  as  they  are  circularly  propagated  nearly  yielding  uniform  electric  field  on 
the  window  surface,  hence  generating  uniform  plasma.  This  structure  produces  argon  plasma  in 
the  pressure  range  of  10-300  mTorr  with  densities  of  the  order  of  101 1  - 1 0 1 2  cm'3,  typical 
microwave  power  being  500  W. 

The  examples  of  structures  we  have  briefly  reviewed  show  that  most  of  them  are  based 
on  surface  wave  excitation  principle,  radiated  or  evanescent  waves  coupling  or  a  combination  of 
them.  In  many  cases,  the  structure  is  a  multimodes  one  and  it  finally  operates  on  a  specific  mode 
depending  on  the  matching  conditions  or  the  resonant  conditions  which  are  used.  In  any  way,  we 
have  seen  that  the  design  of  the  structure  is  really  depending  on  the  wanted  plasma  configuration 
accounting  for  the  aimed  application. 

Nevertheless,  regarding  the  main  characteristics  of  these  structures,  we  have  observed,  in 
all  cases,  that  a  threshold  power  is  required  in  order,  at  least,  to  generate  a  plasma  filling  the 
chamber.  This  threshold  particularly  depends  on  the  power  balance,  say,  the  0/N  parameter 
where  0  is  the  mean  power  for  electron- ion  pair  maintaining  and  N  is  the  neutral  density. 
Concerning  the  plasma  homogeneity,  the  overcoming  difficulty  is  the  general  inhomogeneity  of 


Fig.  1.1 8.  Coaxial-type  open-ended  dielectric  cavity  [79]. 


the  electric  field.  Therefore,  many  innovative  devices  have  been  conceived  to  yield  nearly 
homogeneous  distribution  of  the  electric  field.  The  plasma  homogeneity  is  then  obtained 
decreasing  the  pressure  toward  the  low  pressure  range  provided  the  plasma  diffusion  length  is 
larger  than  the  inhomogeneity  length  of  the  field. 


CHAPTER  II 

PHYSICAL  MODEL  OF  THE  MICROWAVE  DISCHARGE 
ON  A  DIELECTRIC  ANTENNA  SURFACE 

2.1.  Introduction 

The  microwave  discharge  on  a  dielectric  antenna  surface  being  flown  about  can  be  used 
for  boundary  gas  layer  modification,  changing  gas  flow  in  the  antenna  vicinity.  For 
understanding  of  the  physical  processes  taking  place  in  such  kind  of  the  discharge  the 
mathematical  model  of  the  surface  microwave  discharge  should  been  worked  out  and  the 
theoretical  analysis  should  ensure  the  following  results: 

(•)  the  calculation  of  the  dispersion  curves  of  a  surface  wave  taking  into  account 
inhomogeneity  of  a  plasma  layer  (self-consistent  solution  with  the  ionized  non-linearity); 

(••)  the  calculation  of  a  field  penetration  depth  in  plasma  and  thickness  of  a  plasma  layer 
created  by  both  direct  and  reflected  surface  waves,  in  dependence  on  wave  amplitude  and 
frequency,  and  also  gas  pressure  and  airflow  velocity; 

(•••)  the  dependence  of  longitudinal  length  of  the  discharge  supported  by  a  surface  wave 
from  power,  transferred  by  a  microwave; 

(••••)  the  analysis  of  requirements  of  electromagnetic  waves  reflection  from  antenna 
butt-end  at  the  different  shape  of  the  latter; 

(•••••)  the  analysis  of  two  dimensional  field  structure  in  a  plane  of  the  antenna. 

The  geometry  of  installation,  in  which  the  microwave  discharge  is  created,  is  given  in 
Fig.2.1.  A  microwave  energy  forming  plasma  is  supplied  to  one  edge  of  a  dielectric  antenna.  The 
direction  of  the  supersonic  gas  flow  reverses  the  direction  of  the  microwave  spreading.  At  the 
initial  moment  the  plasma  forms  by  strong  microwave  fields  at  the  boundary  between  waveguide 
and  antenna.  Then  the  plasma  starts  to  propagate  along  the  antenna  covering  only  its  part  surface 
at  small  power  (Fig. 2. 2),  or  its  entire  surface  at  high  microwave  power  (Fig. 2. 3).  Electron 
plasma  density  (marked  by  4  in  Fig. 2. 2  and  Fig. 2. 3)  in  a  steady-state  mode  exceeds  critical 
density  that  enables  an  electromagnetic  field  to  propagate  along  plasma  as  a  surface  wave.  The 
electromagnetic  field  penetrates  into  plasma  on  distance  of  the  order  c/coPe  (c  is  the  light 

velocity,  cope  is  the  electron  plasma  frequency),  therefore  plasma  exists  as  a  flat  sheet  near 
antenna.  The  sheet  thickness  depends  on  microwave  wave  power,  energy  transfer  and  diffusion 


Fig.2.1.  Installation  for  creation  of  microwave  discharge  on  antenna  surface. 

1  -  Dielectric  wave-guide  antenna 

2  -  Waveguide  carrying  microwave 


MW 


Fig. 2. 2.  Structure  of  surface  discharge  at  small  power  of  microwave. 

1  -  Dielectric  wave-guide  antenna 

2  -  Waveguide  carrying  microwave 

3  -  Non-  ionized  gas 

4  -  Space,  occupied  by  plasma. 


Fig. 2.3.  Structure  of  surface  discharge  at  high-power  of  microwave. 

1  -  Dielectric  wave-guide  antenna 

2  -  Waveguide  delivering  microwave 

3  -  Weakly  ionized  gas  region 

4  -  Plasma  region 


coefficients  in  plasma.  At  a  large  distance  from  the  antenna  the  electric  field  absents,  therefore 
electrons  intensively  disappear  due  to  attachment.  In  this  area  (marked  by  3  in  a  Fig. 2. 2  and 
Fig. 2. 3)  the  electrons  density  is  small  (usually  plasma  in  this  area  forms  by  photoionization). 

In  surveyed  model  the  gas  flow  in  the  antenna  vicinity  is  considered  to  be  given.  The 
discharge  influence  to  gas  flow  (for  example,  owing  to  a  neutral  heating  by  discharge)  is  not 
taken  into  account. 


2.2  Complete  set  of  equation  for  discharge  description 


2.2.1.  Maxwell  equations 


The  mathematical  model  grounded  on  set  of  Maxwell  equations 

rotE  = - — , 

c  dt 

(2.1) 

-  IdE  4tt  -. 

rotB  =  — —  H - j , 

c  dt  c 

(2.2) 

divE  =  47tep, 

(2.3) 

divB  =  0, 

(2.4) 

and  equations  of  charged  and  neutral  particles  kinetics.  Here  E  and  B  are  electric  and  magnetic 
field  intensity,  j  is  electric  current  density  and  p  is  density  of  a  space  charge.  We  shall  consider 
that  the  frequency  of  microwave,  supporting  plasma,  is  much  greater  then  the  frequency  of 
inelastic  collisions  and  reverse  time  of  energy  exchange  between  electrons  and  neutrals  and 
character  size  of  a  microwave  field  inhomogeneity  is  much  greater  then  thermal  conductivity 
length  kTe  =  ke(M/2m)^2 ,  where  m  and  M  are  electronic  and  ion  masses,  Xc  is  the  mean  electron 

free  path.  In  this  case  we  can  consider  that  electrons  temperature  will  be  local  function  of 
coordinate  and  it  is  quasistationary  function  of  time.  The  plasma  can  be  considered  as 
quasistationary  medium  with  time,  but  without  a  spatial  dispersion,  with  slow  time-dependent 
dielectric  permittivity,  that  means  0  In  (ne)/r)t  ~  V;  «  co.  (V;  is  ionization  frequency).  Thus,  the 
microwave  field  in  plasma  will  be  harmonic: 

(2.5) 


E(r,t)=  E(f,t)exp(-icot), 

B(r,t)  =  B(f,t)exp(-icot). 

Therefore  the  equations  (2. 1-2.4)  can  be  rewritten  as  following  equation: 
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dive(o),r)E  =  0, 
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divB  =  0, 


(2.40 


where  £(©,?)=  1 - 47inc(l)e — _  ,  Vgn  _  coiiision  frequency  between  electrons  and  charged 

mco(co+iven) 


particles,  e>0  -  elementary  electric  charge,  ne(r)  -  density  of  electrons. 

The  similar  approach  is  usual  for  microwave  discharge  calculation  [1-6],  the  large  list  of 
paper  is  given  in  [2],  [24]. 

The  alternate  approach,  grounded  on  a  numerical  integration  of  Maxwell  equations 
during  a  field  period  [5],  essentially  increases  volume  of  calculation  and  is  justified  in  a  case 
when  it  is  necessary  to  take  into  account  non-linear  electromagnetic  field  interactions  (for 
example,  microwave  field  rectification  in  space  charge  layer  on  plasma  boundary).  This  method 
can  be  used  for  discharge  modelling  in  specific  regimes. 


2.2.2.  Kinetics  of  charged  particles  without  negative  ions 


To  determine  of  the  spatial  distribution  of  plasma  density  it  is  necessary  to  solve  the 
kinetics  equations  for  electrons,  ions  and  neutral  excited  particles.  As  it  was  shown  these 
equations  are  essentially  various  in  plasma  (where  characteristic  length  of  inhomogeneity  more 
than  Debye  radius)  and  in  weakly  ionized  gas  (where  the  reverse  condition  takes  place). 

In  first  case  we  have  to  take  into  account: 

1.  Direct  and  step  ionization  of  atoms  (molecules)  by  electron  impact, 

2.  Atoms  and  molecules  exciting  by  electron  impact, 

3.  Ionization  at  heavy  particles  collisions, 

4.  Dissociative  recombination. 

The  time-space  evolution  of  charged  particles  density  is  determined  by  the  balance 


equation  (  ne  =  ni  =  n  ) 

^+div(nCs)-VDaVn  =  F(n.Te). 


(2.7) 


The  equation  (2.7)  describes  plasma  density  distribution  in  all  space  except  for  a  space 
charge  layer  and  Knudsen’s  layer  near  to  antenna  surface.  In  the  formula  (2.7)  the  denotation  for 
an  ambipolar  diffusion  coefficient 


pv  Fe^-t-  +  M-+Pe 

fie +1-4 


(2.8) 


is  used.  The  electron  Ge=neVe  and  ion  G+=n+V+  current  satisfy  the  equations. 


G  =Ge-G+=  0.  (2.9) 

Ge=G+=neCs-DaVne. 

The  relation  (2.9)  assumed,  that  plasma  has  no  electric  current,  that  can  be  created  only 
by  external  electrodes.  In  (2.7)-(2.9)  and  below  the  denotations  Ve,  V+,  Cs  are  velocities  of 
electrons,  ions  and  neutrals,  pe  and  p+  are  mobilities  of  electrons  and  ions,  De  and  D,  are 
coefficients  of  a  free  electron  and  ion  diffusion,  e  is  elementary  electric  charge,  E  is  electric  field 
strength  in  plasma,  F  is  expression  for  sources  of  electron  and  ions  which  are  taking  into  account 
that  they  are  bom  only  by  pairs.  The  items,  containing  Cs,  take  into  account  motion  of  a  gas 
stream  concerning  antenna  resulting  in  to  ions  and  electrons  drift  downstream.  The  standard 
expression  for  F  is 

F  =  -  anen+ ~Pngn+ .  (2.10) 


The  expression  (2.10)  takes  into  consideration  a  direct  ionization  by  electronic  impact  (Vi 
is  frequency  of  ionization),  and  two-  and  three-particles  recombination  (a  and  (1  are 
recombination  coefficients).  Frequency  of  ionization  and  the  recombination  rates  are  functions  of 
temperature  of  electrons.  At  recording  of  equation  (2.9)  we  have  limited  by  slow  processes  with 
typical  frequencies  less  then  ion-neutral  collision  frequency.  It  is  possible  because  the  character 
time  of  electrons  density  change  is  determined  by  ionization  frequency  and  it  is  us  ually  on  some 
order  below. 


2.2.3.  Kinetics  of  charged  particles  with  negative  ions 

The  set  of  equations,  used  in  2.2.1,  satisfactorily  describes  processes  in  steady  discharge 
at  high  electron  density.  For  the  description  of  formation  of  a  transversal  profile  of  electron 
density  and  discharge  propagation  along  the  antenna  the  kinetics  of  electrons  at  small  density, 
which  should  take  into  account  attachment-detachment  reaction,  will  play  essential  role. 


The  majority  of  reaction  including  negative  ions  can  be  described  by  the  scheme  1 


e  +  A  +  B  <->  A  +  B  . 

Constants  of  a  forward  reaction  (attachment)  and  reverse  reactions  (detachment)  we 
designate  ka  (dimension  -  cm3/s)  and  kd  (cm6/s).  Except  for  attachment  and  detachment  we  take 
into  account  also  recombination  positive  and  negatively  ionized  atoms  (reaction  constant  (L 
(cm6/s))  Density  of  negatively  charged  ions  we  designate  n_  and  their  current  density  -  G_ 

In  the  given  paragraph  we  use  also  v,  (cm3/s)  and  (3  (cm6/s)  as  notations  for  ionization 
constants  by  electronic  impact  and  three-body  recombination 

e  +  A  <-^e  +  e  +  A+ . 

Thus  initial  set  of  equations  will  accept  a  view: 

3nP 


3t 

3nH 

Ik 

3n_ 

“aT 


-  +  divGe  =  (vjN  -(3nen+)ne  -  (kaNne  -kdn_)N  , 

-  +  divG+  =  (vjN  -(3nen+)ne  -(3_nen_, 

-+  divG_  =  (kaNne  -kdn_)N  -(3_nen_ . 


(2.11) 

(2.12) 

(2.13) 


Electron  Ge=  neVe  and  positive  and  negative  G+=  n+V+  and  G_=  n_V_  ion  flows  fit  to 


equations 


Ge=  neCS  -BeneE  “  DeVne  > 

(2.14) 

G+=n+Cs  +  [t+n+E  -  D+V  n+ , 

(2.15) 

G_=  n_Cs  -  |l_n_E  -  D_Vn_ , 

(2.16) 

j  =  eG  =  e  (Oj-Gg  -G_ ) , 

(2.17) 

divE  =  -47te(ne  +  n_  -  n+ ) , 

(2.18) 

rotE  =  0  .  (2.19) 

In  this  system  we  take  into  account  attachment  and  detachment  of  electrons  (constant  of 

reaction  ka  and  kd)  and  also  recombination  of  positive  and  negatively  ionized  atoms. 


1  As  an  example  it  is  possible  to  record  reaction  [42] 

e  +  Ot  +  O2  ^  O2  +  O 
e  +  Ot  +  Nt  >  Oo  +  N 


e  +  O+Oo  > O  +0 


c  +  O  +  Nt  <-^0  +N 


In  the  region  of  small  electron  densities  (rDe  »  L),  the  equations  (2. 1 1)— (2. 1 8)  yield: 


On 


dt 

On 


On 


dt 


=  0, 

(2.20) 

-Vcp, 

(2.21) 

■  +  div(neCs  +peneVcp 

)-DeAne  =  (ViN-pnen+)ne  -(kaNne  -kdn_)N, 

(2.22) 

-  +  div  (n+Cs  -  (i+n+Vcf 

>)“  D+ An+  =  (vjN  -  pnen+  )ne  -  p_n+n_ , 

(2.23) 

-  +  div  (n_Cs  +  |i_n_V  cf 

>)-  D_An_  =  (kaNne  -kdn_)N  -(3_n+n_ . 

(2.24) 

In  a  reverse  limiting  case  ( rDe  «  L ,  ne  +  n_  =  n+ ,  div  G  =  0 )  we  have 

div  {((|ie  +  p+  )ne  +  (|!_  +  (l+  >i_  )Vcp}-  (De  -  D+  )Ane  -  (D_  -  D+  )An_  =  0 , 

E  =  -Vcp, 

j  =  eG  =  e(-  ((|ie  +  (i+  )ne  +  (|le  +  (l_  )n_  )Vcp  +  (De  -  D+  )Vne  +  (D_  -  D+  )Vn_) , 


5n-  +  div(neCs  +(ieneVcp)-DeAne  =  (v^ -[3nen+)ne  -(kaNne  -kdn_)N  , 
+  div  (n_Cs  +  |l_n_Vcp)-  D_An_  =  (kaNne  - kdn_ )N  -  (3_ (ne  +  n_)n_ . 


dt 

dn 


dt 


(2.25) 

(2.26) 

(2.27) 

(2.28) 
(2.29) 


On  the  contrary  in  this  case  both  ambipolar  field  and  field  supporting  a  direct  current 
appear  not  potential  (contain  a  vortex  component),  though  a  general  electrical  field  is  potential. 
Therefore  equation  (2.25)  includes  expressions,  keeping  the  second  derivative  of  density  of 
charged  particles.2 


The  alternate  approach,  founded  on  the  introduction  of  a  current  vector-potential  does  not  result 
in  simplification  of  equations  because  the  vortex  electric  current  can  exist  in  plasma  without 
external  electrodes,  that  follows  from  equation  (2.30).  In  this  case  set  of  equations  becomes 


rot 


rotA  +  (De  -  D+)Vne  +  (D_  -  D+  )Vn 
(|le+(l+)ne +  (ji_+ (i+)n_ 


=  0, 


(2.30) 


^  +  div  (6e  )  =  (VjN-  (3ne  (ne  +  n_  ))ne  -  (kaNne  -  kdn_)N  ,  (2.31) 

^-^  +  div(G_)=  (kaNne  -  kdn_)N  -  P_(ne  +n_)n_,  (2.32) 

dt 

g  _  rotA  -  (Dc  -  D+)Vne  -  (D_  -  D+  )Vn_ 

(M-e  +  M-+  )ne  +  (ft—  +  M-+  )n_ 


(2.33) 


The  simplification  of  a  set  of  equations  is  possible,  if  in  plasma  the  electron  conduction 
predominates  (|ie  +  |a+)ne  »  (|i_  +p+)n_,  or,  that  practically  is  equivalent,  ne  »  (m/M)^2n_. 
In  this  case  at  the  description  of  an  ionic  flow,  neglecting  an  ionic  diffusion  and  conductivity,  we 
receive 


div  {(fiene  )V\|/}=  0 ,  (2.37) 

E  =  -Vy-— Vlnne,  (2.38) 

M'e 


at 


+  div 


neCs  +  (ien  eV^/  ■ 


(ji+n+  -(l_n_)De 


Vn 


!Gne 


=  (vjN  -  (3ne  (ne  +  n_  ))ne  -  (kaNne  -  kdn_  )N 


3n+ 

~dt 

5n_ 

~dt 


-  +  div 


n+Cs  -|l+n+V\)/- 


fVVPe 

l^ene 


Vm 


-  +  div 


V 

'  H-n-De 


=  (v|N  -  (3nen+  )ne  -p_n+n_ , 


(2.39) 


(2.40) 


) 

\ 


n_Cs  +  |l_n_V\|/  +  ^ 


Vn„ 


M-e^e 


(kaNne-kdn_)N-(3_(ne  +n_)n_ .  (2.41) 

v  r '^c  J 

In  a  reverse  limiting  case  (pc  +  p+  )ne  «  (|i_  +  p+)n_  ambipolar  condition  bind  motion  of 

positive  and  negative  ions,  at  the  same  time  electrons  moves  as  free  particles,  as  the  reference 
value  of  an  ambipolar  potential  will  be  equal  to  k(T+  +  T_)/e  .  The  set  of  equations  will  accept  a 
view: 

div  (n(p_  +  p.+  )grad\|/)  =  0  ,  (2.42) 


r  =nr _ \wcotA 

0-^e  M-+  )nc  0-1 —  T  M-+  )n- 

(|lc  D+  +  |l+De  )nc  +  ({!_  +  q+  )Dcn_ 

(|le+(i+)ne  +  ((i_  +  (i+)n_ 
(l_n_rotA 


Vo.  + 


_ M-e^eO^-  P+) _ 

((le +lt+)ne +((!_  + (t+ )n_ 


Vn_ 


G_=  n_Cc  - 


([le  +  H+  )ne  +  (|l_  +  (l+  )n_ 
(lt_D+  +  [i+D- )n_  +  (de  + 11+  )D_n 


(lie  +ll+)ne  +  (jll_  +|Ll+)n_ 


—  Vn  + 


(l_n_(De-D+) 


(lte+li+)ne  +  (li— +  lt+)n_ 


Vn„ 


G+=  (ne  +  n_)Cs  +  ,  Mn.+n)ro.A 

(llc  +  [1+ )ne  +  ()l_  +  (l+  )n_ 

(lteD+  +(1+De)nc  +  (lt+Dc  +lt_D+)n_ 

((le  +  lt+)ne  +  (|1_ +lt+)n_ 

(lteD+  +(i+D-)ne  +  (|a_D+  +  li+D_)n_  ^ 
(ltc  +(t+)ne  +  (|i_  +  q+)n_ 


(2.34) 


(2.35) 


(2.36) 


(2.43) 


E  =  — V\|/-  — — — Vhn+  , 

F-+F+ 


3n, 


at 


■  +  div 


n+Cs  -(l+n+V\)/- 


|i+D_  +(i_D4 


•  Vn , 


v  H-+H+  j 

=  (vjN  -  (3ne  (ne  +  n_  ))ne  -  (3_n+n_ 


3n_ 

~aT 


+  div 


n_Cs  +|x_ilVv|/- 


(i+D_  +(i_DH 
(i_  +  (i+ 


Vn_ 


=  (kaNne-kdn_)N-(3_(ne  +  n_)n_ 


(2.44) 


(2.45) 


^  +  div  (neCs  -  DeVne  )=  (vjN  - pne  (ne  +  n_))ne  -  (kaNne  -kdn_)N  .  (2.46) 


The  description  of  electron  detachment  kinetics  is  principle  for  the  description  of 
discharge  processes.  Detachment  is  connected  to  interaction  with  exited  atoms  and  molecules, 
and  also  with  heating  of  negatively  ionized  atoms  in  an  electrical  field  (in  a  corona  discharge). 
The  fullest  models  of  a  kinetics  used  for  the  description  of  streamers  propagation  include  up  to 
150  [7-10]  of  chemical  reactions.  Features  of  our  system  is  the  availability  of  two  areas  of 
discharge.  The  first  area  is  the  area  of  weak  ionization.  In  this  area  the  excited  atom  density  is 

small,  and  thus  detachment  is  unessential  kd^  =  0.  The  second  area  is  area  of  high  density  of 
electrons.  In  this  area  the  electron  and  excited  atoms  concentrations  are  great.  Approximately  it 
is  possible  to  consider,  that  the  electrons,  excited  atoms  and  negatively  ionized  atoms  are  in 
equilibrium  in  this  area. 


neN  _  §N  f 

n-  g- 


2;tkTe  6 


m 


exp 


kT.. 


(2.47) 


where  £_  is  the  bond  energy  of  electron  in  negative  ion. 

At  application  of  a  principle  of  detail  equilibrium  the  ratio  (2.47)  allows  to  calculate  a 
detachment  constant 


,(2)  =  v  =  k  Sn 

vd  Ka  Ka 

n_  g- 


27tkTe  V2 
m 


exp 


(  £_  ^ 


kT(. 


(2.48) 


At  small  electron  concentrations  the  excited  atoms  are  formed  by  electronic  impact, 
therefore  N*~kne,  and  detachment  probability  kd^  =  kne/N.  With  growth  of  an  electron 
density  the  number  of  electrons  borne  due  to  detachment  should  increase  from 

dN^/dt  ~  knen_N  up  to  dN^/dt  ~  kd2V_N  .  In  the  given  report  we  model  this  transition  by 


the  introducing  of  a  phenomenological  function  <b(ne)  such,  that  lim  <b(ne)  =  0, 

ne— >0 

lim  <E>(ne)  =  l.  Detachment  rate  we  shall  model  by  expression 

ne 

dN/dt  =  dN(l)/dt  (l  -  4>(ne ))  +  dN(2)/dt  <f>(ne ) . 

The  complete  set  of  chemical  reactions  in  air  can  be  found  in  papers  [7-14].  In  the  initial 
equations  set  we  can  neglect  a  space  charge  in  area  of  weak  ionization,  therefore  these  equations 
can  be  reduced  to  the  form 


> 

-e 

ii 

o 

(2.49) 

5te  +divGe  =  Fe(ne,n+,n_,Te), 

(2.50) 

+  +divG+  =  F+(ne,n+,n_,Te), 
dt 

(2.51) 

+  divG_  =  F_(ne,n+,n_,Te), 
ot 

(2.52) 

Ge=  neCs  +  [teneVcp  -  DeVne , 

(2.53) 

G+=  n+Cs  -p+n+Vcp-D+Vn+, 

(2.54) 

G_=  n_Cs  +p_n_Vcp-D_Vn_, 

(2.55) 

E  =  -Vcp . 

(2.56) 

The  subsequent  perturbation  on  small  parameter  r\  1  =(rDe/L)  1  can  be  received  with  the 
equation 

Acp(1)  =  47te(ne -n+ +  n_).  (2.57) 

The  requirements  of  particle  fluxes  through  boundary  and  ambipolar  field  potential 
continuity  are  executed  on  boundary  between  weakly  ionized  gas  and  plasma.  On  plasma  and 
antenna  boundary  the  space  charge  layer  takes  place.  The  boundary  conditions  can  be  noted  in 
the  form  of  the  equality  of  electrons  and  ions  currents  on  wall 

jen  J-n  —  jin  ’  (2.58) 

and  Bohm  criterion  [15,  16] 

ln  =  nivs-  (2-59) 

The  frequencies  of  chemical  processes  are  determined  by  electron  and  neutral 
temperatures  defined  by  energy  balance  equations. 


2.2.4.  Electron  energy  balance  and  calculation  of  electron  temperature 


As  follows  from  equation  (2.9),  in  a  limiting  case  of  strongly  ionised  gas  (L»rDe, 


ne  =  n+  =  n )  the  electrical  field  exists  in  two  forms:  ambipolar  (longitudinal)  field 

E  =  -  °e  ~  D+  V  ln(n)  =  ^kTe/e-p+kT+/e -in(n)^ 

He  +  H+  He  +  H+ 

and  field  that  supported  electric  current,  flowing  past  on  plasma,  (transversal  field) 
G 


E  = 


-  =  -V\|/. 


(He+H+)n 

In  double-fluid  hydrodynamics  approximation  the  electron  energy  balance  equation 
reduced  to  [39,  40] 


3  dnkT 
2  dt 


+  div 


3 


— nkTV 

V2  7 


+  nkTdiv  V  +  divq  =  Q+  -  Q_ . 


(2.60) 


The  terms  in  the  left-hand  part  of  the  equation  (2.60)  take  into  account  the  action  of 
electron  pressure.  In  a  right-hand  one  Q+  is  heat  energy  transmitted  to  electrons  by  an  electrical 
field  (action  of  the  conforming  friction  force),  Q_  is  thermal  energy  transmitted  from  electrons 
to  the  neutrals  by  collision. 

Power  of  an  electrical  field  transmitted  to  electrons,  can  be  counted  as  work  of  an 
electrical  field  in  coordinate  system,  bound  with  driving  neutral  gas 

Q+e=-e((Ge-nCs)E).  (2.61) 

Substituting  the  expressions  for  electric  field  and  current  (2.9)  in  (2.60)  we  shall  receive: 


Q+e  — 


HeG2  Da(De-D  +  )(vnf  (GVn)peDe  -  2peD+  -p+Dc 


(2.62) 


n(|ic+H+)2  (He+H+)  n  n  (He  +  H+  )2 

It  is  more  convenient  to  use  expression  which  is  taking  into  account  a  smallness  of  ion 
temperature  in  comparison  to  electron  one. 


HeG2  H+D2  (vn)2  (GVn)  [leDe 


£+e 


(2.62a) 


n  (He  +  H+)2  (ne  +  H+)2  n  n  (He  +  H+  )2 
The  first  term  in  equation  (2.62)  responds  for  electrons  heating  by  electric  current, 
flowing  past  on  plasma,  the  second  one  accounts  energy  expended  for  ambipolar  field  creation 
by  electrons,  the  last  term  is  cross  addends  which  are  taking  into  account  electron  current 
through  an  ambipolar  barrier.  In  microwave  field  it  is  necessary  to  add  into  equation  (2.62)  the 


n  e2v  E2 

power  of  high  frequency  field  transferred  to  electrons  (including  both  collisional  — e(  2en 


(o32 +vjn)’ 


2mlco  +v 

and  collisionless  absorption  of  energy)  (see  for  example  [1]  -  [5]). 

The  similar  expression  can  be  recorded  for  energy  transferred  from  electric  field  to  ions 


£++  ' 


|1+G2  Da(De-D+)(vn)2  (GVn)|ieD+ +  2|i+De -|i+D, 

/■>  ‘  /  \  t  •>  r\ 


(2.63) 


n(|ie  +  |i+)2  (f G  +  B+)  n  n  (lte+l02 

The  conforming  simplified  expression  which  is  taking  into  account  the  smallness  of 
electron  temperature  looks  like 


Q 


++ 


[l+G2  [i+D2  (Vnf  (GVn)  2[i+De 

,  r\  '  ✓  r\  , 


(2.63?) 


n(|ae  +  |i+)2  (|ie+!t+)2  n  n  (|ae+M-+)2 

As  well  as  it  was  necessary  to  expect,  the  energy  spent  by  electrons  for  creation  of  an 
ambipolar  barrier,  is  transmitted  to  heating  of  ions,  so  that  as  a  whole  ambipolar  field  does  not 
work,  it  only  reallocate  energy  between  components  of  plasma.  The  effect  of  ions  acceleration 


1  The  indicated  result  can  be  illustrated  by  following  calculation.  Let  the  electron  is  in  an 
electrical  field  E  and  let  it  experiences  only  isotropizing  collisions  t1,t2,---,tn  in  points  with 
coordinates  q ,  r2 ,  •  •  • ,  rn .  Then 

=7i-i  +^j-i(tj-ti-i)+|^-(tj-tj-i)2’ 

m(V/-Vn). 

2  v  J_1  /VJ  J-1'  2m 

Substituting  in  last  expression  Vj_  j  from  first  and  summarizing  on  all  interferences  we 


e(vj_,EXtj-tj_1)+A-(tj-tHk 


receive 


mfe-Vp) 


e((rn~%)E). 


Thus,  as  well  as  it  was  necessary  to  expect,  the  change  of  a  thermal  energy  of  an  electron 
is  equal  to  change  of  its  potential  energy 
3  3kT^ 


—  n - 

2  dt 


where  the  current  j  includes  all  possible  causes  of  electron  movement  both  the  current  caused  by 
an  electrical  field  (conduction  current)  and  current  conditioned  by  inhomogeneous  distribution  of 
electrons  in  space  and  their  Brownian  motion  (diffusive  current).  Allowing  follow  up  elastic 
energy  losses  of  electrons  in  interferences,  the  equation  of  electron  energy  balance  can  be 
reduced  to  the  form  (thermodiffusion  and  heat  conduction  can  be  taken  easily  into  account) 

( (  2  3  3 

ne 


3  dkT 

2 n  at 


mv 


-E  +eDeVn 
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E 

J  J 


--Kvk(Te-Tg). 


by  an  ambipolar  field  is  well-known  in  low  pressure  discharge,  restricted  by  walls,  where  this 
ions  energy  is  transferred  to  a  wall.  In  freely  localized  discharges  this  energy  transmitted  to 
heavy  neutral  particles  due  to  ion  neutral  collisions.  The  estimations  demonstrate,  that  this 
energy  can  play  a  noticeable  role  in  heating  heave  plasma  components4 5. 

Similar  calculation  for  not  quasi-neutral  area,  where  T|  =  rDe/L»l  for  electrons  and 
similar  ratio  for  ions  results  in 

Q+e  =fienee(vcp)r  -De(VneVcp),  (2.64) 

Q++  =|i+n+e(vcp)2  -D+(vn+Vcp).  (2.65) 

The  electrical  fields  in  this  area  can  be  great,  and  density  of  charged  particles  is  small, 
therefore  we  can  neglected  the  second  term  in  equations  (2.64),  (2.65). 

The  calculation  of  power  transmission  from  electrons  in  a  other  component  should  take 
into  account  the  ionization  losses,  excitation  of  electronic,  vibratory  and  rotary  levels  and  elastic 
losses.  All  losses  we  can  take  into  consideration  by  introducing  an  effective  energy  transfer 
coefficient  K  [30].  In  this  case  expression  for  Qe_  has  the  form 


Qe— =  KVenne-(kTe-kTg). 


(2.66) 


Thus,  the  electron  temperature  in  a  quasi-stationary  mode  can  be  calculated  under  expression: 
2(Qe+-Qe_-WP) 


T  =T  +  - 

j.e  i 


3kKVenne 


(2.67) 


where  Wp  is  energy  lost  due  to  ambipolar  field  creation,  spotted  by  the  left-hand  part  of  an 


equation  (2.60)/ 


This  yields  the  expressions  (2.62)  and  (2.62?).  The  physical  sense  of  diffusive  cooling  of 
electrons  is,  that  the  majority  of  electrons  located  in  a  point  r  were  born  and  acquire  energy  in 
the  field  of  large  integrated  electron  densities,  therefore  having  come  exactly  to  r ,  they  have 
expended  a  part  of  the  thermal  energy  on  overcoming  of  a  potential  hill. 

4  It  follows  from  mass  equality  of  ions  and  neutrals.  The  heating  of  the  neutrals  at  transition  of 
ionization  front  can  be  estimated  approximately  under  the  expression 

AT  =  ApTe , 

where  A  -  factor  about  units  dependent  on  a  spacing  distribution  of  chemical  reactions  at  the 
front  of  ionization,  p  -  degree  of  ionization. 

5  At  record  of  equation  (2.67)  we  have  determined  the  temperature  T  as  3/2  squares  of  mean 
energy  in  coordinate  system,  bound  with  driving  electrons,  thus  temperature  is  determined  by  an 
isotropic  part  of  a  cumulative  distribution  function  in  this  coordinate  system.  The  ionization 
frequency  is  determined  through  a  cumulative  distribution  function  in  coordinate  system  bound 
with  neutral  gas,  therefore  we  should  characterize  mean  electron  energy  in  this  coordinate 


2.3.  Surface  and  dielectric  waveguide  waves  in  homogeneous  bounded  plasma 


The  dispersion  characteristics  of  surface  waves  were  esteemed  in  number  of  papers. 
Among  the  24  works  of  Russian  scientists  it  is  necessary  to  mark  the  monographies  [17-23].  The 
fullest  review  of  the  articles,  concerning  surface  wave  discharge  structure  investigation  keeping 
more  of  160  links,  was  published  by  Zhelyazkov,  and  Atanassov  [24].  The  fullest  subsequent 
reviews  of  recent  works  are  contained  in  stuffs  of  conferences  [25,  26]. 

For  surface  wave  dispersion  curve  calculation  in  a  considered  case  it  is  convenient  to  use 
the  approach  offered  in  [27].  Let's  consider,  that  the  surface  wave  is  propagated  in  z-direction, 
and  the  plasma  is  inhomogeneous  only  in  a  direction,  perpendicular  to  its  boundary  (this 
direction  we  shall  designate  X).  As  it  is  known,  the  field  of  a  surface  wave  has  three  components 
Ex  ,  Ez  and  By  .  Let's  search  for  the  solution  by  the  way: 


r  > 

Ex 

Ex(x) 

Ez 

>  =  < 

Ez(x) 

By 

By(x) 

exp  (- i(cot  -  ihz )) . 


(2.68) 


In  this  case  Maxwell  equations  are  resulted  in 
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(2.70) 


(2.71) 


In  range  of  parameters,  interesting  us,  we  can  leave  out  a  spatial  dispersion.  Therefore 


e(x)  =  l  — 


47tnc‘ 


(2.72) 


mco((0+iven) 

The  boundary  conditions  for  an  equation  (2.69)-(2.72)  include  conditions  of  radiation 
non  reflection  from  infinity  and  conditions  of  electrical  and  magnetic  field  continuity  on  plasma 


system  (conforming  quasi-temperature  we  shall  designate  as  Te  =  Te  +  m  /?> )  And  this  mean 

energy  will  determine  frequencies  of  chemical  processes,  bound  with  interferences  of  electrons 
and  neutrals.  It  is  determined  by  an  equation 
2Qe+ 


T  =T  + 

e 


3kKVenne 


(3.10) 


Also  is  in  the  full  consent  with  an  equation  of  an  energy  balance  recorded  for  one  electron. 


and  dielectric  antenna  boundary.  The  field  inside  the  antenna  can  be  shown  by  the  way 

•exp  (-/(cot  -ihz)). 
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(h  cl  COE  D  )ch  {\j{h~  —co  £d/ c2)(.r  —  x0 )) 
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(2.73) 


where  Xo  is  x- coordinate  of  a  central  plane  of  the  antenna. 

If  the  electron  oscillation  frequency  in  an  electric  field  in  a  direction  perpendicular  to 
the  border  exceeds  the  typical  size  of  a  non-uniformity  of  a  space  charge  layer  (electronic  radius 
of  Debye)  the  process  of  UHF  field  rectification  in  a  space  charge  layer  and  the  process  of  the 
harmonic  generations,  which  are  capable  to  update  process  of  absorption  of  energy  in  plasma, 
are  essential. 

The  electromagnetic  fields  in  such  system  can  be  separated  into  three  groups  (Fig.  2.4). 

1.  The  waves  of  an  continuum  spectrum,  radiated  by  system  at  lr<co2£p/c2  (if  £p<0  this 
waves  damp  along  the  antenna  and  exist  only  near  to  sharp  structural  inhomogeneities  -  at 
transferring  point  from  a  waveguide  to  the  antenna  and  near  to  a  point  of  plasma  slab  abruption). 

2.  The  waves  of  a  dielectric  wave  guide  (their  field  is  periodic  on  X  direction  inside 
waveguide  and  damps  exponentially  outside).  The  expression  for  the  dispersion  characteristics  of 
these  waves  has  no  analytical  form  for  arbitrary  plasma  density,  however  at  defined  values  of  £p 
the  wave  propagation  coefficient  easily  evaluate  (m=0,l,2...). 
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(The  wave  for  which  hn2<k02£P  does  not  exist). 

3.  The  surface  wave  existing  at  £p<-£d.  The  surface  wave  concentrates  near  to  plasma 
and  antenna  boundary  and  exponentially  damps  in  both  sides  from  it. 
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The  penetration  depth  of  a  surface  wave  and  waves  of  a  dielectric  waveguide  in  plasma 
determines  by  relations 


L=Im  k 
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kx2=ko2£p-h2. 
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Fig.  2.4.  The  waves  in  system  a  plasma  -  dielectric  waveguide  in  absence  of  collisions. 

1  -  The  surface  wave  existing  in  an  area  £p<-£d; 

2  -  dielectric  waveguide  wave  existing  in  an  area  £p>£d; 

3  -  boundary  infinite  spectrum,  4  -  k02£o; 

5  -  calculation  under  the  equation  (2.7)  (eD  =  5  ,  antenna  thickness  L=8  mm, 
microwave  length  X=2.4  cm,  nc  =  mco2/47te2  ). 
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Fig.  2.5.  Inverse  penetration  depth  x 1  of  a  surface  wave 
as  a  function  of  electron  density  at  v/co=0. 
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The  calculations  results  for  surface  wave  and  dielectric  waveguide  waves  for  geometry 
used  in  experiment  (L=l  cm)  are  given  in  Fig.  2.4. 

The  dependence  of  return  depth  of  penetration  of  a  surface  wave  in  plasma 
corresponding  to  a  dispersion  curve  in  Fig.  2.4  is  submitted  on  Fig.  2.5.  In  area  e;,  — >  — eD  it  can 
be  described  by  equation 


_i  co 
x  =  — 


+  £r 


An  utilized  earlier  estimation 


-1  ®pe 
X  =  — — 

C 

is  take  place  at  ne  — >  °° . 

Maximum  UHF  field  penetration  depth  in  plasma 

..-l.^pe  2£D 
c  £d  —  1 

is  watched  at  £p  =  -2ed  or,  that  is  the  same,  ne  =  nee(l  +  2ed  ). 

At  nc  — >  nec  (ed  + 1)  the  depth  of  penetration  of  a  field  in  plasma  decreases  up  to  zero 
point  (at  absence  of  collisions)  and  at  definite  density  of  plasma  becomes  equal  or  more  than 
thickness  of  a  diffusive  layer.  Therefore  at  ne  =nee(l+2£D)  the  considered  model  requires  in 
refinement.  It  is  possible  to  expect,  that  in  this  area  the  electromagnetic  wave  will  be  absorbed 
much  more  strongly,  than  follows  from  our  model. 

The  electrons  collisions  result  in  attenuation  both  surface  and  waveguide  modes,  at  that 
the  surface  wave  exists  in  area  Re  h  >k0  Re  £d,  and  it  turns  in  strongly  attenuated  wave  with  a 
reverse  dispersion  at  kfRe  £p<Re  h2<k02Re  £d-  The  results  of  calculation  of  the  surface  wave 
dispersion  characteristics  at  the  account  of  collisions  are  given  in  Fig.  2.6. 


2.4.  Surface  wave  discharge  classification 


A  microwave  discharge  on  a  surface  of  a  dielectric  wave-guide  antenna  can  be 

■j  o 

supported  in  pressure  range  10"  -10  torr.  In  general  case  such  wide  pressure  range  results  in 
necessity  of  the  account  of  all  possible  kinetic  processes: 

1.  Direct  and  step  ionization  by  electron  impact  and  ionization  at  collision  of  excited 
molecules. 


Fig.  2.6.  The  surface  wave  dispersion  characteristics  in  system 
of  a  plasma  -  dielectric  waveguide  at  the  account  of  collisions  v/O)=0,01. 

9  9 

Surface  wave:  1  -  Re  h  ,  2  -  Imh“.  Continuum  spectrum  boundary  -  3. 


2.  Attachment  and  detachment. 

3.  Electron  -  ion  and  ion  -  ion  recombination. 

4.  Excitation  of  the  oscillatory  and  rotary  freedom  degree  of  molecules. 

5.  Ambipolar  diffusion. 

6.  Particles  and  energies  balance  can  be  both  local  and  non-local  depending  on  external 
conditions. 

7.  At  low  pressures  it  is  necessary  to  include  a  spatial  dispersion  in  electrodynamics’  model. 

8.  The  discharge  can  be  transparent  and  opaque  (at  high  pressures)  for  self-radiation. 

9.  The  plasma  non-uniformity  along  an  electrical  field  direction  results  in  additional  energy 
absorption  in  a  resonance  region,  where  the  real  part  of  complex  permittivity  of  plasma  is 
close  to  zero  point.  There  is  a  considerable  change  of  an  electron  distribution  function 
due  to  additional  energy  absorption  in  a  resonant  layer. 

10.  Sharp  changing  of  the  electric  field  near  plasma  boundary  can  change  the  spatial 
distribution  of  the  plasma  density,  and  caused  the  generation  of  the  electric  field 
harmonic  due  to  striction  plasma  nonlinearity. 

1 1 .  The  role  of  the  gas  heating  can  be  different  depending  on  a  part  of  energy  transmitted  to 
the  heating  of  neutrals  and  on  a  microwave  pulse  duration’s. 

The  complete  set  of  the  constant  rates  of  chemical  reactions  for  air  plasma  is  given  in  [7- 

10]. 

We  shall  use  notifications  below:  e  and  m  -  electron  charge  and  mass;  c  -  light  velocity; 
ne  -  charged  particles  density;  ven  -  electron  collision  frequency;  co  -  microwave  field  frequency; 
A.en  -  mean  free  path  of  electron-neutral  collision;  CQpe=(47tnee2/me)1/2  -  plasma  frequency; 

nec  =  meco2/47te2  -  critical  density;  kQ  =  co/c  ;  E  and  B  -  electric  and  magnetic  fields;  e(co)  - 

dielectric  permittivity;  in  antenna  region  e(co)  =  £D;  in  plasma  e(co)  =  £P  =  1  - (Qpe/ ot(ot  +  iven ); 
Da  -  ambipolar  diffusion  coefficient;  v,  -  ionization  frequency;  a  -  recombination  coefficient; 

cy  -  thermal  capacity;  x  -  electron  thermal  conductivity  (heat  transfer)  coefficient;  Te  and  Tg  - 
electron  and  gas  temperatures.  For  estimation  of  field  penetration  length  into  plasma  we  will  be 
used  the  value  c/ cope  as  for  surface  microwave  discharges  cope  >  co . 

As  the  position  of  one  of  the  discharge  boundaries  does  not  fixed  by  a  wall,  the  usual 
analysis  based  on  the  solution  of  the  ambipolar  diffusion  equations  and  direct  -  current  analogy 
can  not  be  used.  In  the  report  we  shall  give  approximate  model  for  stationary  microwave 


discharge  valid  at  middle  gas  pressure  (p=l-100  torr)  and  on  its  base  the  classification  of  the 
main  processes  taking  place  in  plasma  of  surface  microwave  discharge  will  be  given. 

If  the  condition  c70Jpe<A.cn  is  satisfied  the  skin  layer  on  plasma  boundary  is  collisionless, 
therefore  processes  of  energy  transfer  from  microwave  field  to  electrons  should  be  considered  in 
kinetic  approach.  Otherwise  the  hydrodynamic  equations  are  possible  to  use.  The  border  of  the 
kinetic  and  hydrodynamic  description  on  the  diagram  (ne-p)  -  (charged  particles  density  - 


pressure  of  the  neutrals)  is  described  by  a  curve 


,  2^2  P  ’ 
47te  Ael 


(2.75) 


where  Ael  is  the  mean  free  path  of  electrons  at  pressure  p=l  torr. 

The  energies  balance  is  local  if  the  characteristic  size  of  absorption  of  microwave 
energy  c/tOpe  for  typical  case  exceeds  the  length  of  thermal  conductivity  )iTe=^en/(S)1/2,  where  8  is 
a  mean  part  of  energy,  which  electrons  loss  in  collisions  with  molecules  (at  elastic  collisions 
S=2m/M,  that  yields  a  lower  limit  of  possible  values).  The  equation  of  border  of  local  and  non¬ 


local  modes  looks  like 
c2mS 


c  mo  2 
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(2.76) 


Particles  balance  is  local  at  condition  arie  >D(cope/c)“  for  the  local  energies  balance  case 
and  anc>D8/(Aen)2  in  the  opposite  case.  The  boundary  position  between  these  areas  on  the 
diagram  (ne  -  p)  results  in  equations  (Di  is  an  ambipolar  diffusion  coefficient  at  pressure  1  torr) 


D,  47te 
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D,  8 

ne  =  9  p. 

«  4i 


(2.77) 


(2.78) 


In  the  case  of  a  three-particle  recombination  predominance  (third  particle  is  an  electron, 
P  is  a  recombination  coefficient)  the  equations  (3.3)  and  (3.4)  become 
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Thus  the  electrons  density  and  temperature  will  be  connected  by  Sakha  equation. 


The  conducted  above  reasoning  are  fair  for  plasma  by  transparent  for  own  radiation, 

otherwise  the  energy  absorbed  by  electrons  leaves  the  discharge  volume  only  by  heat  conduction 

(including  radiant  heat  conduction),  and  the  cooling  of  discharge  descends  only  from  a  surface. 

The  Rosseland’s  path  of  a  photon  can  be  estimated  under  the  expression 

AT7/2 
L  =  — , 


where  the  value  of  a  constant  A  depends  on  a  kind  of  gas  and  on  a  photons  scattering  mechanism 
[28],  [29].  As  the  value  ^en/VS  is  the  typical  size  of  discharge  at  local  energies  balance  and 
c/(Qpe  is  the  typical  one  at  non-local  energies  balance  the  conditions  of  the  discharge  opacity 
look  like 
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(2.80) 


If  the  conditions  (2.79)  or  (2.80)  are  fulfilled  the  discharge  will  be  equilibrium.  The  discharge 
transition  from  the  non-equilibrium  form  to  equilibrium  usually  occurs  in  time  as  development  of 
the  overheated  instabilities  [30]. 

The  neutral  gas  heating  in  discharge  area  results  in  two  effects. 

1.  There  is  a  hydrodynamic  gas  movement  leading  to  equalising  of  a  gas  pressure  in  the 
discharge  with  the  pressure  of  the  environmental  medium.  The  typical  time  x  of  this  process  is 
f  7/2 
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for  the  systems  with  local  energy  balance  and 

c  \[d 

x  ~  x - P 

^elQs 

for  the  system  with  non-local  one. 


(2.82) 


2.  At  large  times  the  gas  heating  leads  to  slow  gas  convection  due  to  Archimedes  force. 
This  class  of  phenomena  we  do  not  consider  in  this  work. 

The  areas  with  different  character  of  discharge  kinetic  processes  depending  on  electron 
density  and  the  gas  pressure  are  given  on  Fig.  2.7. 

At  all  conditional  character  of  similar  division  due  to  a  proximity  of  the  estimations  of  a 


Fig.  2.7.  Different  discharge  kinetics  areas: 

I  is  area  of  kinetic  viewing, 

II  is  local  particles  balance  and  non-local  energy  balance, 

III  is  local  particles  and  energy  balance, 

IV  is  non-local  particles  balance  and  non-local  energy  balance, 

V  is  non-local  particles  balance  and  local  energies  balance, 

VI  is  the  area  of  the  discharge  opaque  for  optical  radiation. 


Fig.  2.8.  Different  discharge  kinetics  areas: 

I  is  area  of  kinetic  viewing, 

II  is  local  particles  balance  and  non-local  energy  balance, 

III  is  local  particles  and  energy  balance, 

IV  is  non-local  particles  balance  and  non-local  energy  balance, 

V  is  non-local  particles  balance  and  local  energies  balance. 

VI  is  region  of  discharge  opaque  for  optical  radiation. 

VII  is  the  boundary  of  100%  ionisation. 

VIII  is  the  critical  density  (for  2.45  cm  microwave). 

IX  is  the  boundary  of  field  absorption  in  resonance  point  and 
collisionless  field  absorption  on  plasma  boundary, 

X  is  the  boundary  of  elementary  homogeneous  plasma  model. 
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field  penetration  depth  and  dependence  of  the  chemical  reactions  constants  rates  and  part  of 
energy,  loss  electrons,  on  temperature,  this  analysis  allows  to  make  a  deduction  about  possible 
methods  of  the  plasma  description.  According  to  the  carried  out  analysis  in  large  range  of 
requirements  it  is  possible  to  consider  that  the  particles  balance  is  local,  but  at  the  same  time,  it  is 
necessary  to  take  into  account  the  electronic  thermal  conductivity. 

The  analysis  of  the  discharge  processes  we  shall  begin  from  most  simple  (and  at  the 
same  time  most  often  meeting)  particular  case  of  discharge  with  local  balance  of  the  particles 
and  non-local  balance  of  energies.  On  the  first  stage  we  also  shall  not  take  into  account  the 
hydrodynamic  movement  of  gas  due  to  its  heating. 

Let's  consider  that  the  longitudinal  plasma  inhomogeneity  is  small  and  the  distribution 
of  electromagnetic  wave  along  the  plasma  -  antenna  surface  it  is  possible  to  describe  in  a 
geometric  optics  approximation  [7,8].  It  is  usual  condition  for  discharges  supported  by  a  surface 
wave,  generated  near  to  one  of  the  plasma  borders.  As  the  area  of  the  microwave  field 
localization  in  plasma  is  much  less  than  thermal  conductivity  length  it  is  possible  to  account  that 
the  energy  absorption  distribution  in  a  direction,  perpendicular  to  a  boundary,  is  described  as 

W±=W±0exp(-x^Aj.  (2.83) 

Considering  that  the  part  of  energy  expended  on  ionization  does  not  depend  on  energy  (it 
means,  that  electron  density  is  proportional  to  absorbed  energy,  that  executed  only 
approximately),  it  is  possible  to  receive  following  expression  for  dependence  of  electrons  density 
on  cross  coordinate 

ne(x) =  neo  exp  (-  x^/A,en ),  (2.84) 

where  neo=Vi(W_L0)/oc.  Thus  the  characteristic  size  of  plasma  exceeds  the  depth  of  electric  field 
penetration  in  plasma  therefore  from  an  electrodynamic  point  of  view  it  is  possible  to  consider 
plasma  to  be  uniform  with  electron  density  equals  to  its  boundary  value.  Strictly  speaking,  on 
plasma  and  dielectric  antenna  boundary  the  electron  density  is  close  to  zero,  therefore  the 
diffusive  layer  is  formed  in  this  area. 

2.5.  Simplified  set  of  equations 

2.5.1.  Electron  heat  transfer  equation 

The  homogeneous  plasma  model  can  be  use  for  calculation  of  surface  wave  dispersion 
characteristics,  if  the  reference  size  of  a  diffusive  layer  is  mach  shorter  then  the  surface  wave 
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penetration  depth  in  plasma.  The  matching  demonstrates  that  the  plasma  non-uniformity  should 
essentially  change  surface  wave  spectra,  at  first,  in  the  area  of  small  electron  densities 

n  c2  v 

|eD|  <  £p  <  2|ed|,  and  secondly,  in  the  area  of  high  densities  —  >  — — — .  In  intermediate  area 

nc  co  Da 

2 

2eD  <  —  <  — — —  the  elementary  model  must  work,  except,  maybe,  the  account  of  absorption 
nc  or  Da 

in  a  resonant  point  [30],  [31].  The  boundary  of  this  area  is  showed  in  a  Fig.  2.8  by  curves  X. 
Thus,  the  model  of  homogeneous  plasma  works  presumptively  only  in  the  field  of  high-pressure 
discharge,  at  pressure  below  than  30  torr,  one  has  to  use  the  numerical  models. 

At  a  discharge  numerical  modelling  one  can  outline  some  possible  problems: 

1.  Simulation  of  longitudinal  discharge  structure  (along  a  direction  of  a  wave  propagation). 

2.  Simulation  of  discharge  parameters  distribution  on  the  whole  antenna  surface. 

3.  Simulation  of  stationary  discharge  structure  in  a  direction,  perpendicular  to  antenna 
boundary. 

4.  Simulation  temporary  evolution  of  transversal  discharge  structure. 

The  independent  consideration  of  "longitudinal"  and  "transversal"  processes  can  be  done 
because  the  longitudinal  discharge  size,  as  a  rule,  is  much  greater  then  its  transversal  size  and 
thus  stationary  longitudinal  structure  is  reached  much  earlier  than  the  longitudinal  one.  At  the 
same  time  the  discharge  time  evolution  modelling  should  include  the  hydrodynamic  motion  of 
neutral  gas,  which  one  can  include  its  convective  motion.  This  motion  depends  on  a  plenty  of  the 
factors,  including  arrangement  of  the  antenna  in  a  horizontal  or  vertical  plane.  Therefore 
simulation  of  process  of  installation  of  fixed  frame  of  discharge  in  many  respects  has  conditional 
nature. 

As  follows  from  the  previous  analysis  for  discharge  mathematical  modelling  under  our 
condition  (pressure  of  the  neutrals  10-1000  torr,  microwave  wavelength  2.4  cm)  it  is  necessary  to 
take  into  account  following  processes: 

1.  Electron  energy  transfer  described  by  an  equation  of  electrons  energy  balance 

|nk  ^  =  V%VTe  +  Q(x)-  Wl (x)n ,  (2.85) 

2  dt 

where  X  is  electronic  heat  conduction  coefficient;  Q(x)  is  power  transmitted  from  a  surface  wave 
to  electrons,  which  include  both  collisional  (ohmic)  and  collisionless  mechanisms;  wi  is  energy 
lost  by  an  electron  in  interferences,  including  elastic  losses,  loss  on  excitation  and  ionisation 
[30].  Let's  mark,  that  we  consider  discharge  to  be  transparent  for  radiation,  behind  exception, 
may  be,  resonance  emission. 


CLC\ 


2.  Excited  atoms  and  molecules  density  transfer,  including  resonant  condition, 
described  by  an  integral  Biberman-Holstein  equation,  and  the  metastable  atoms,  described  by  a 
usual  diffusion  equation.  It  is  necessary  to  account  excitation  transfer  separately  from  electrons 
energy  transfer  (though  the  some  energy  levels  can  be  in  equilibrium  with  electrons)  because 
region,  where  the  electronic  heat  conduction  is  possible,  is  limited  to  a  potential  barrier  built  by 
an  ambipolar  field,  while  area,  where  excitation  transfer  in  neutral  gas  is  possible,  is  not  limited. 

3.  The  electrons  transfer  in  space  is  described  by  an  ambipolar  diffusion  equation 
(2.7)  taking  into  account  the  ionization  from  exited  states 

Distribution  of  an  electromagnetic  field  in  space  we  shall  describe  by  equations  (2.2)  - 
(2.4).  Energy  absorption  in  a  resonant  point  or  collisionless  absorption  of  an  electromagnetic 
wave  on  plasma  boundary  we  shall  allow  by  the  introducing  of  the  conforming  surface 
resistance. 


2.5.2.  Electron  density  distribution  near  antenna  surface 

In  §2.6  the  different  types  of  UHF  discharge  in  neighbourhood  of  a  dielectric  antenna 
(see  Fig.  2.3)  are  analyzed  and  the  approximated  qualitative  model  of  the  discharge  supported  by 
a  surface  wave  is  constructed.  This  model  used  the  assumption  of  local  charged  particles  balance 
and  a  non-local  energy  balance  in  plasma.  This  model  corresponds  to  the  most  widespread 
discharge  maintenance  conditions  (Fig.  2.7). 

The  surface  waves  characteristics  used  in  model,  leave  out  of  the  corrections,  bound  with 
a  plasma  non-uniformity  near  to  a  dielectric  antenna  and  electromagnetic  field  absorption  in  a 
resonance  region,  where  Re  8  =  0 ,  and  also  possible  non-linear  processes  in  a  space  charge  layer 
between  plasma  and  antenna. 

At  calculation  of  electron  density  distribution  near  plasma  boundary  we  will  use  the 
approach  of  ambipolar  diffusion.  Thus  the  electron  density  distribution  is  subject  to  an  equation 

— -VDaVn  =  Vin-an2 -pn3.  (2.85) 

8t 

The  notations  are  entered  in  (2.85):  n  is  density  of  charged  particles,  Da  is  coefficient  of 
an  ambipolar  diffusion.  The  expression  (2.85)  takes  into  consideration  a  direct  ionization  by 
electron  impact  (Vi,  is  frequency  of  ionization;  a  and  (1  are  two  and  three  body  recombination 
coefficients).  Frequency  of  ionization  and  the  frequencies  of  a  recombination  are  functions  of 
temperature  of  electrons.  As  we  consider  processes  in  a  diffusive  layer,  which  has  size  of  the 


order  •x/Da/vj  .  It  size  is  much  less  than  heat  conduction  length  ^Te=^en/^-  we  shall 
consider  electron  temperature  and,  therefore,  constant  of  chemical  reactions  not  dependent  from 
co-ordinates. 

This  equation  should  be  resolved  with  boundary  conditions  on  infinity  and  on  plasma 
boundary  with  a  dielectric  antenna.  The  former  is  n  — >  nQ  at  x  — »  ,  where  nD  is  the  root  of 
equation 

Vi-an^-pn^O.  (2.86) 


The  elementary  version  of  boundary  conditions  on  a  wall  n  =  0  at  x=0  was  suggested  by 
Shottky.  However  we  must  determine  the  plasma  density  on  the  boundary  with  space  charge 
layer  (CPL).  We  shall  use  a  condition  following  from  B ohm's  criterion 

nVs=D,|h  (2.87). 

At  a  solving  the  equation  (2.87)  it  is  most  convenient  to  consider  that  the  maximum  value 
of  an  electron  density  is  r\„  entering  a  new  variable  T|  =  n/nG  .  Thus  in  a  stationary  approach  the 


equation  (2.85)  is  reduced  to 

2 

— r  =  r|a(l  -  r|)  +  j3  (l  -  rf ) 
dx“ 

where  the  notations 


a  =  on0/Da,  p  =  «n2/Da 


are  entered.  The  equation  (2.88)  can  be  integrated 


dr| 

dx 


a  B  R  4  2a  3 
— ri  + — ri 

3  2  2  3 


(2.88) 


(2.89) 


(2.90) 


The  reference  depth  of  a  diffusive  layer  can  be  evaluated  as 

x  =  v/Da/Vi=(a  +  p)'1/2.  (2.91) 


In  Schottky-type  boundary  condition  the  second  integrating  leads  to  simple  relation 
between  electron  density  and  coordinate 


■  l 

=1- 


dr| 


la  B  R  4  2a  3 
— +  — +  — ri  + — rr-  a  + 
3  2  2  3 


(2.92) 


Bom's  criterion-type  boundary  condition  (2.87)  yields 


/CO 


hs  —  +  ^-+  ^-n4  +  —  rj3  -(a  +  p)r|2 
V  3  2  2  3 


(2.93) 


The  value  T|v  is  determined  from  a  condition  (2.87)  together  with  an  equation  (2.90). 


Thus,  we  have 


vsns  =  Da^+|+|ns  +  ^i-(a+p)ns  • 

Thus,  the  electron  density  on  boundary  is  determined  as  the  root  of  equation 

a  P  p  4  2a  3  (~  X  V§  A 

—  +  —  +— r\s+ — tic  —  a  +  B+— ns  =  0. 

3  2  2  ^  3  1  D; 

v  a  7 

Considering  T|s  «  1 ,  we  shall  receive  approximately 


(2.94) 


(2.95) 


a  p 
ric  =  — +  — 

s  A  3  2 


oc  +  p  +  ■ 


(2.96) 


After  substitution  of  values  a  and  P  from  a  ratio  (2.89)  we  shall  receive: 


an„  pn“  /  o  2  Vc  ann  pn~  l  „  ?  )  ~~ 

~T  +  /  ano  +  PnS  +  tt-  =  ,  ~T+^r  /lan0 +  Pno  +  vin)'  (2.97) 


3  2 


where  Vin  -  ion  neutral  collision  frequency. 


3  2 


Last  ratio  demonstrates,  that  the  plasma  resonance  area  [31,  32]  exists  in  plasma  if 


conditions 


no  >  nC  = ' 


(2.98) 


mot 

Tlsno  <  nc  =  - — 2 
47te 


(2.99) 


is  fulfilled,  where  ot  -frequency  of  microwave  field  supporting  plasma.  The  ratio  (2.99)  results  in 
following  connection  of  electron  density  and  pressure  p 


+  / [ano  +  Pno  +  Viop)n0  • 


(2.100) 


If  condition  (2.100)  is  fulfilled,  additional  surface  wave  energy  absorbs  in  the  plasma 
resonance,  which  results  in  surface  wave  dispersion  curve  change  (increase  of  absorption),  and, 
therefore,  change  the  relations  between  electron  density  and  wavelength,  that  leads  to  more 
strong  dependence  of  electron  density  from  longitudinal  coordinate.  It  gives  the  possibility  to 


support  discharge  at  smaller  microwave  powers.  The  boundary,  which  separates  area  where  the 
resonance  area  is  present  at  plasma,  showed  by  curves  VIII  and  IX  in  Fig.  2.8.  Thus  the 
microwave  length  was  considered  equal  to  2.4  cm.  Above  the  curve  IX  the  plasma  resonance  is 
inside  of  a  spatial  charge  layer,  therefore  the  collisionless  absorption  of  a  surface  wave  can  be 
taken  into  account  in  this  case  as  absorption  on  plasma  boundary,  similar  to  absorption  of  a  HF 
field  in  a  capacitive  HF  discharge  [15].  The  discharge  with  density  is  less  then  critical  one 


nc  =  mco2/47te2  (below  by  curve  VIII)  as  a  rule,  will  not  be  realized1. 

Except  for  curves  VIII  and  IX  Fig.  2.8  contains  the  straight  line  VII  representing  limit  of 
100  %  ionization.  Strictly  speaking,  all  reduced  estimations  are  fair  much  below  by  curve  VII. 
Besides  Fig.  2.8  contains  a  position  of  boundaries  of  discharge  with  different  types  of  kinetics 
(curves  I- VI),  which  position  and  physical  sense  were  discussed  in  §2.2. 

At  low  electronic  density,  when  the  double  recombination  dominates,  it  follows  from 
(4.18)  that  the  influence  of  a  diffusive  layer  on  the  dispersion  characteristics  of  an 
electromagnetic  field  is  essential  at 
,2 

(2.101) 
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that  results  in  a  condition  coincidence  with  a  condition  (2.76),  however  having  sense  not  only  on 
border  of  III  and  V  areas  but  also  inside  II  area.  In  this  case  the  qualitative  analysis  of  discharge 
characteristics  reduced  in  the  given  section  is  inapplicable  and  the  numerical  solution  of  a 
problem  is  necessary. 


2.6.  Discharge  characteristics 


The  equations  set  (2.1/)-(2.4/)  and  (2.7)-(2.9)  were  solved  in  several  cases. 

(•)  Hie  dispersion  curves  of  a  surface  wave  taking  into  account  inhomogeneity  of  a 
plasma  layer  (self-consistent  solution  with  the  ionized  non-linearity)  were  calculated.  In  this 
case  set  of  equations  were  solved  on  segment  (0<y<L)  as  eigenvalues  problem  for  a  nonlinear 
wave  propagated  along  the  antenna.  Solution  of  a  problem  is  searched  as 
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'ky) 

ex 


p(ihz)' 
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(2.102) 


1  Usually,  to  create  discharge  with  density  below  critical  it  is  necessary  to  use  discharge  in  a 
resonator,  when  even  the  little  density  change  results  in  sharp  decreasing  of  a  microwave  field 
due  to  resonance  frequency  shift. 


Boundary  conditions  for  a  field  are  requirement  of  symmetry  at  y  =  0  and  the 
requirement  not  reflections  of  a  wave  in  a  point  y  =  L  (which  with  adequate  accuracy  by  virtue 

of  a  field  shielding  can  be  exchanged  zero).  The  eigenvalue  of  the  task  is  propagation  constant  h. 

The  similar  problem  for  a  wave,  supporting  the  discharge  in  the  waveguide,  was  solved  in 
article  [33],  where  the  discharge  tube  was  introduced  in  a  rectangular  waveguide.  The  analogous 
calculation  for  a  surface  wave  discharge  was  fulfilled  by  Moisan  and  Zakrzhevsky.  In  distinction 
from  the  indicated  works,  in  our  case,  one  of  plasma  boundaries  is  free,  i.e.  its  position  should  be 
found  from  calculation.  The  similar  problems  were  considered  in  a  Boev’s  cycle  of  papers  [34- 
37]  and  paper  [38],  where  the  plasma  dielectric  permittivity  was  considered  as  local  nonlinear 
function  of  a  high  frequency  field.  In  these  papers  the  field  solutions,  looking  as  wave  channels 
propagating  along  plasma  waveguides  (created  by  power  microwave  itself),  were  found. 
Unfortunately,  direct  comparison  of  results  of  the  articles  [34-37]  with  experiment  is  impossible 
as  the  authors  have  not  taken  into  account  the  electrons  temperature  and  density  space  relaxation 
processes  (i.e.  diffusion  and  thermal  conductivity). 

In  a  nonlinear  task  on  eigenvalues  (opposite  linear  one)  the  value  of  propagation  constant 
h  appears  to  be  the  function  of  field  amplitude.  Calculation  of  dependence  h  =  h(E0),  as  well  as 
the  calculation  of  the  profile  of  the  transversal  distributions  of  plasma  density  n  =  n(E0,y)  and 
fields  E  =  E(E0,y)  ,  is  a  main  problem. 

(••)  The  surface  wave  attenuation  constant  usually  has  a  peak  near  to  resonance  point 
maf 

n  = - T(l +8^),  therefore  at  low  powers  the  surface  wave  is  well  absorbed  and  the  reflection 

4m~ 

from  back  plasma  boundary  is  insignificant  (Fig.  2.8).  In  this  case  it  is  possible  approximately  to 
consider  that  the  plasma  is  supported  by  a  surface  wave,  running  along  antenna  from  microwave 
source. 

At  large  powers  the  plasma  covers  whole  antenna.  The  antenna  butt-end  became  the 
short-circuited  plunger  and  it  is  possible  to  expect  the  essential  reflection  of  a  surface  wave. 
Such  reflection  can  change  a  cross  profile  of  created  plasma,  and  it  is  essential  to  affect  on 
efficiency  of  energy  transfer  from  magnetron  to  plasma  and  magnetron  operation  regime.  The 
energy  deposition  to  a  neutral  gas  also  becomes  periodic  function  of  longitudinal  coordinate. 
Under  these  conditions  the  problems  (•••)  and  (••••)  is  planed  to  consider  analytically  and  then 
to  include  them  in  numerical  discharge  model.  One  of  the  problems  is  the  analysis  of  a  surface 
wave  reflection  from  butt-end  of  the  different  shape  antenna  with  the  purpose  to  achieve 
minimum  wave  reflection  coefficient. 


The  propagated  surface  wave  energy,  absorbed  per  unit  plasma  length,  can  be  expressed 
through  an  energy  flux  density  along  the  antenna 

Wabs  =2WImh(ne)  (2.103) 

that  allows  by  virtue  of  a  requirement  (2.102)  to  calculate  power  necessary  for  maintaining  of  the 
surface  discharge  plasma  with  density  ne : 


W  -  neWl 
2Imh(ne) 


(2.104) 


where  wi  is  the  power  expended  for  creation  of  one  free  electron.  The  calculated  dependencies 
(2.104)  in  relative  units  for  the  different  values  of  ratio  v/co  of  collision  frequency  v  to 
microwave  field  frequency  to  is  represented  in  Fig.  2.9. 

The  joint  solution  of  the  equation  (2.104)  and  Maxwell  equations  allows  to  receive 
approximate  dependence  of  plasma  parameters  on  the  propagated  surface  microwave  power 


Distribution  of  surface  wave  power  along  a  plasma  slab  one  can  find  from  the  equation: 
dW(z) 


— —  =  -2  Im(h(ne(W(z)))  )w(z) 
dz 

The  equation  (2.105)  is  equivalent  to  the  equation  for  electron  density: 

dne(z)  _  2Im(h(ne(z))) 

dZ  ln(W(ne(z))) 

dm 


(2.105) 


(2.106) 


where  the  dependence  W(ne(z))  is  shown  in  Fig.  2.9.  We  neglect  wave  reflection  from  back  butt- 
end  of  plasma.  It  is  obviously  possible  if  the  wave  does  not  reach  the  butt-end  of  the  antenna. 

The  calculated  dependencies  of  plasma  density  on  coordinate  Z  along  the  antenna  are 
given  in  Fig.  2.10.  In  summary  on  Fig.  2.11-2.14  the  results  of  calculation  of  distribution  of  a 
very  high  frequency  field  in  neighbourhood  vicinity  of  the  antenna  at  different  integrated 
electron  densities  in  plasma  ne/nec  are  submitted.  The  existence  of  a  diffusive  layer  on  border  of 
plasma  is  not  taken  into  account.  As  well  as  it  was  necessary  to  expect,  at  ne  — >  °°  field 
allocation  comes  nearer  to  allocation  in  a  metal  wave  guide. 
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Fig.  2.9.  Electrons  density  to  electrons  critical  density  ratio  rie/riec  as  a  function  of  power  W  (in 
relative  units)  of  a  surface  wave  propagated  along  the  antenna  at  different  collision  frequencies 


v/co:  1-0.01;  2-0.03;  3-0.1;  4-0.3;  5-1. 


Fig.  2.10.  Electrons  density  to  electrons  critical  density  ratio  ne/nec  as  a  function  of  coordinate  Z 
at  different  frequencies  of  collisions  v/co:  1  -  0.01;  2  -  0.03;  3  -  0.1;  4  -  0.3;  5  -  1.0. 


ne/nec=2,5 

Fig.  2.11.  A  space  distribution  of  electrical  (Ey,  -  curve  (1)  and  Ez  -  curve  (2)) 
and  magnetic  (By  -  curve  (3))  fields  in  a  surface  wave. 


ne/nec=5  (greatest  depth  of  penetration) 

Fig.  2.12.  A  space  distribution  of  electrical  (Ey,  -  curve  (1)  and  Ez  -  curve  (2)) 
and  magnetic  (By  -  curve  (3))  fields  in  a  surface  wave. 
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Fig.  2.13.  A  space  distribution  of  electrical  (Ey,  -  curve  (1)  and  Ez  -  curve  (2)) 
and  magnetic  (By  -  curve  (3))  fields  in  a  surface  wave. 
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Fig.  2.14.  A  space  distribution  of  electrical  (Ey,  -  curve  (1)  and  Ez  -  curve  (2)) 
and  magnetic  (By  -  curve  (3))  fields  in  a  surface  wave. 
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CHAPTER  III 


EXPERIMENTAL  INSTALLATIONS 


3.1.  Experimental  set-up 


The  experimental  set  up  is  based  on  a  stainless  steel  cylindrical  gas/vacuum  chamber 
three  meters  long  with  about  a  meter  in  diameter  (Fig. 3.1).  Supersonic  airflow  can  be  organized 
in  it.  Twenty  six  windows  with  different  diameters  (10-50  cm)  are  situated  on  its  lateral  surface, 
that  has  allowed  to  visualize  the  discharge  processes  and  to  provide  a  vacuum-tight  transport  of 
electric  current  and  electric  signals.  The  gas/vacuum  chamber  consists  of  two  sections  that  can 
form  a  vacuum-tight  junction  with  help  of  a  lever-operated  gate.  At  installation  and  adjustment 
of  diagnostic  and  experimental  equipment  inside  the  chamber  its  sections  can  be  easily 
separated.  A  vacuum  pump  has  allowed  to  achieve  a  pressure  level  of  1  torr  for  10-15  minutes  of 
pumping. 

The  supersonic  airflow  was  formed  at  filling  the  chamber  with  the  atmospheric  air 
through  specially  profiled  converging-diverging  nozzles  with  Mach  number  M=2.  A  detailed 
calculation  of  the  nozzle  shape  has  been  carried  out.  The  nozzle  convergent  part  was  made  close 
to  the  ideal  form  providing  a  plane  M=1  interface.  The  walls  near  the  critical  cross  section  had  a 
sharp  edge.  The  angle  of  the  wall  of  the  divergent  part  near  the  critical  cross  section  was 
calculated  from  the  theoretical  formulae  with  use  of  the  gas  dynamic  functions.  The  whole 
profile  of  the  divergent  part  was  calculated  to  obtain  the  output  parallel  supersonic  airflow.  This 
uniform  M=Mo  airflow  region  in  case  of  negligible  pressure  in  the  vacuum  chamber  is  restricted 
by  a  surface  with  a  shape  of  two  Mach  angle  cones  based  on  the  nozzle  exit  circle;  upstream 
M<Mo,  downstream  M>Mo.  The  nozzles  were  made  of  a  dielectric  (caprolactone).  A  cylindrical 
nozzle  had  an  outlet  diameter  3  cm  and  a  diameter  15.4  mm  of  critical  cross  section.  The 
characteristic  period  of  time  of  regular  nozzle  operation  was  about  2  second. 

The  nozzles  were  screwed  on  a  dielectric  tube  5  cm  in  diameter.  On  its  opposite  end  a 
electric  hydraulic  valve  with  05  cm  channel  was  fixed.  All  this  construction  was  installed  on  a 
flange  of  a  window  of  the  chamber  so  that  to  position  the  supersonic  airflow  along  the  chamber’s 
diameter.  The  nozzle  end  was  situated  on  the  level  of  big  transparent  windows  on  the  lateral 
surface  of  the  chamber.  The  electric  hydraulic  valve  operated  with  a  power  supply  unit  which 
was  synchronized  with  the  discharge  power  supply  units. 


Fig. 3.1.  General  scheme  of  the  experimental  set  up 


The  experimental  set  up  was  equipped  with  different  types  of  diagnostics. 

A  Schlieren  device  with  015  cm  outlet  lens  was  used  for  visualization  of  the  free 
supersonic  airflow  and  of  the  discharges  in  it.  An  electric  incandescent  lamp  with  radiating  body 
area  about  1  mm”  was  used  as  a  light  source  for  adjustment  of  the  device.  At  registration  of  the 
airflow  gas  dynamical  structure  a  pulsed  (x=5  (is)  xenon  lamp  was  applied. 

A  spectrograph  STE-1  and  monochromator  MDR-3  were  used  for  registration  of  plasma 
radiation  spectra.  Either  a  highly  sensitive  photo  film,  or  a  CCD  image  sensor  was  fixed  in  the 
focal  plane  of  the  spectral  devices.  For  an  experimental  investigation  of  a  surface  microwave 
discharge  in  supersonic  flow  it  is  necessary  to  use  contactless  methods  and  the  optical 
diagnostics  of  plasma  fully  satisfy  this  requirement.  The  vibrational  temperature  Tv  was 
measured  by  the  relative  intensity  of  the  nitrogen  and  CN  bands.  The  gas  temperature  Tg  was 
determined  by  the  distribution  of  the  intensity  of  rotational  lines  of  molecular  bands  of  a  second 
positive  system  of  nitrogen. 

Double  ray  pulsed  oscilloscopes  with  pass  band  10  MHz  and  sensitivity  1  mV/cm  were 
used  for  registration  of  the  discharge  voltage  characteristics  and  the  electric  photomultiplier  tube 
signals. 

For  obtaining  of  the  threshold  characteristics  of  a  surface  microwave  discharge  the 
dependencies  of  the  minimum  input  power,  at  which  the  discharge  on  an  external  surface  of  a 
dielectric  body  starts  to  be  formed,  on  different  values  of  air  pressure  in  a  vacuum  chamber  and 
on  different  durations  of  microwave  pulses  were  measured.  In  experiments  the  moment  of 
formation  of  the  discharge  was  registered  visually  or  on  appearance  of  a  signal  on  the  screen  of  a 
double -beam  oscilloscope  from  the  collimated  photoelectric  multiplying  tube,  tuned  on  area  of 
the  antenna  at  an  edge  of  a  waveguide.  Thus,  on  the  second  beam  of  an  oscilloscope  the  signal 
from  a  microwave  detector  was  applied.  The  amplitude  of  this  signal  was  proportional  to  a 
pulsed  microwave  power. 

The  common  view  of  a  surface  microwave  discharge  was  registered  on  a  film  and  a  video 
camera  in  two  projections  (side  view  and  top  view).  It  has  allowed  to  fix  the  longitudinal  size  of 
a  surface  discharge  and  to  measure  longitudinal  velocity  of  its  distribution  at  different  air 
pressures,  pulse  durations  and  microwave  powers. 

For  measurement  of  gas  and  vibration  temperatures  under  condition  of  the  surface 
microwave  discharge  in  supersonic  airflow,  a  diagnostic  set  up  has  been  assembled.  Radiation  of 
the  discharge  that  is  formed  in  supersonic  airflow,  is  projected  on  input  slots  of  two  spectral 
devices  through  the  windows  of  chamber,  the  mirrors,  and  the  focusing  quartz  lenses.  The  lens 
diameter  and  its  focal  length  were  chosen  so  that  to  provide  projection  of  the  discharge  image 


with  reduction  k=10,  and  use  of  all  the  light-gathering  power  of  the  devices.  A  spectrograph 
STE-1  and  monochromator  MDR-3  were  used,  their  inverse  linear  dispersions  being  0.35  and 
1.3  nm/mm  accordingly  for  the  spectral  region  of  3000-6000  A0.  The  measuring  system 
(including  the  spectrograph  STE-1,  monochromator  MDR-3  and  the  photo  film)  was  calibrated 
by  registration  of  standard  iron  arc  spectrum  (to  get  a  wavelength  scale)  through  a  9- step  light 
reducer  (to  reveal  the  linear  region  of  film  darkening)  on  the  same  photo  film  together  with  the 
tested  spectrum.  The  photo  films  were  analysed  with  use  of  a  microphotometer  IFO-451. 

The  electronic  analyzer  of  optical  spectra  on  base  of  a  personal  computer  registers  plasma 
radiation  in  visible,  ultraviolet  and  infrared  spectral  ranges.  The  system  is  designed  on  base  of  an 
IBM-compatible  personal  computer.  The  sensor  of  plasma  radiation  on  base  of  a  CCD  image 
sensor  forms  a  video  signal  directly  proportional  to  the  radiation  intensity  in  the  spectral  band 
3000-9000  A0.  The  video  signal  is  coded  by  an  analog-to-digital  converter  situated  on  an 
Interface  card  and  is  sent  to  the  PC  memory  as  an  array  of  data.  The  discharge  spectra  can  be 
visualized  on  the  monitor  or  printed  out  with  use  of  program  codes.  The  system  can  operate  in 
two  regimes:  synchronous  and  asynchronous.  At  synchronous  operation  the  pulses  of  discharge 
are  synchronized  with  the  CCD  cycles  of  exposition  and  pickup.  This  regime  makes  it  possible 
to  research  pulsed  discharges  or  to  integrate  the  information  on  pulsed-periodical  discharges.  The 
asynchronous  operation  yields  spectra  that  are  averaged  during  all  the  time  of  exposition. 

The  video  signal  formed  by  the  CCD  is  passed  to  a  broad-band  differential  amplifier  for 
compensation  of  a  constant  component  and  amplification  to  the  level  sufficient  for  further  coding 
by  the  analog-to-digital  converter.  Drivers  of  digital  signals  are  necessary  for  formation  of  CCD 
controlling  pulses. 

The  interface  card  contains  a  12-bit  analog-to-digital  converter,  a  quartz  generator  and  a 
scheme  of  generation  of  synchronizing  pulses  for  controlling  of  CCD,  a  universal  parallel 
interface  of  control  of  executive  devices  (valves,  pumps  etc.)  and  an  interface  logic  for  fitting  the 
card  with  the  ISA  bus. 

3.2.  Set-up  for  creation  of  surface  microwave  discharge 

The  system  of  creation  of  the  surface  microwave  discharge  in  a  boundary  layer  near  to 
the  bodies  streamlined  supersonic  airflow  was  designed  and  manufactured.  The  modernization  of 
the  microwave  generator  used  for  creation  of  the  surface  microwave  discharge  inside  a  boundary 
layer,  existing  near  bodies  of  different  configuration  streamlined  supersonic  airflow,  is  carried 
out.  The  experimental  set-up  is  completed  with  the  system  of  creation  of  supersonic  airflow;  the 


special  system  of  input  of  a  microwave  energy  into  discharge  chamber  and  transformations  of 
the  microwave  energy  into  the  energy  of  a  surface  wave;  the  system  of  recording  of  the 
parameters  of  microwave  radiation  and  supersonic  airflow,  the  diagnostic  instrumentation  for 
measuring  of  plasma  parameters  of  the  discharge  and  for  recording  of  the  photo  and  video 
information  about  common  view  of  the  surface  microwave  discharge  in  supersonic  airflow. 

The  pulsed  magnetron  generator  of  a  centimetric  wave  frequency  was  used  as  a  source  of 
a  microwave  energy.  The  magnetron  generator  had  the  following  characteristics:  the  wavelength 
A=2,4  cm;  the  pulsed  microwave  power  W<600  kW;  the  pulse  duration  t=1-F300  ps;  the  pulse 
period-to-pulse  duration  ratio  Q=1000. 

A  pulsed  modulator  with  the  partial  discharge  of  capacity,  developed  in  our  laboratory, 
was  used  for  a  feeding  of  magnetron.  The  work  of  the  pulse  modulator  is  based  on  a  principle  of 
rather  long  accumulation  of  energy  during  a  pause  between  pulses  and  short-term  its  feedback  to 
microwave  generator  during  a  pulse.  The  microwave  energy  introduced  into  the  vacuum 
discharge  chamber  by  the  waveguide  system  of  rectangular  section  9,5x19  mm  through  the 
directional  coupler.  All  waveguide  system  was  hermetic.  It  was  filled  by  SF6  at  pressures  up  to 
5  atmospheres  for  avoidance  of  an  electrical  breakdown.  The  discharge  was  formed  in  a 
cylindrical  chamber.  The  vacuum  system  allowed  to  investigate  the  surface  microwave  discharge 
at  air  pressures  p=  1, 0-760  torr.  The  waveguide  entered  into  the  discharge  chamber  ended  the 
special  antenna,  on  which  the  surface  microwave  discharge  in  supersonic  airflow  was  created. 
The  discharge  was  formed  on  an  external  surface  of  a  dielectric  body  of  rectangular  section 
S=lx2  cm  and  length  L=15  cm.  The  direction  of  supersonic  flow  was  opposite  to  the  direction  of 
the  surface  microwave  discharge  spreading. 

The  output  signal  from  the  microwave  detector  was  registered  on  an  oscillograph.  The 
form  of  this  signal  was  close  to  rectangular.  The  same  signal  acted  on  an  input  of  the  pulse 
digital  voltmeter,  under  which  indications  the  amplitude  of  a  microwave  pulse  was  determined. 
By  meanings  of  this  amplitude  the  microwave  power  supplied  to  the  discharge  chamber  was 
measured.  For  it,  the  calibration  of  the  voltmeter  with  the  help  of  the  power  calorimeter  was 
previously  made.  During  calibration  the  power  calorimeter  was  connected  to  an  output  of  the 
basic  shoulder  of  the  directional  branching-off. 

The  new  set-up  for  realization  of  the  investigations  of  the  surface  microwave  discharge  in 
large  range  of  gas  pressures  ?=10~3-103  torr  is  designed  and  assembled.  The  set-up  consists  of  the 
vacuum  booster  pump  permitting  to  execute  the  preliminary  pumping-out  of  the  discharge 
chamber  up  to  pressures  -10"1  torr;  the  vacuum  turbomolecular  pump  for  pumping-out  of  the 
discharge  chamber  up  to  pressures  ~10"5  torr;  the  glass  discharge  chamber  with  35  cm  in 


diameter  and  50  cm  in  length.  The  chamber  has  a  metal  flange  for  connection  with  the  pumping 
and  inlet  gas  system.  The  set-up  has  the  system  of  measurement  of  gas  pressure,  the  system  of 
the  supply  of  a  microwave  energy  to  dielectric  bodies,  the  system  of  input  into  the  discharge 
chamber  of  the  dipole  antenna  and  electrical  probe  permitting  to  measure  of  distributions  of  the 
electric  field  and  electron  concentration  in  zone  of  existence  of  plasma  of  the  surface  discharge. 

3.3.  Set-up  for  creation  of  surface  transversal  pulsed-periodical  discharge 

The  supply  unit  for  obtaining  the  pulsed-periodic  discharge  on  a  surface  of  a  plate  placed 
at  supersonic  airflow  is  designed  and  created.  The  parameters  of  the  supply  unit  allow  to  create 
the  surface  discharge  in  supersonic  airflow  at  two  modes: 

1.  A  mode  of  single  pulsed  discharge  at  voltage  U=3-25  kV,  current  i=0,l-40  A,  pulse 
duration  x=l  (is  -  1  ms,  maximal  pulsed  power  WP<1  MW. 

2.  A  mode  of  pulsed-periodic  discharge  at  voltage  U=3-25  kV,  current  i=0,l-10  A,  pulse 
duration  x=  1-300  (is,  pulse  repetition  frequency  f=l-100Hz,  mean  power 
consumption  about  1  kW. 

The  system  of  creation  of  a  pulsed-periodic  discharge  on  a  surface  of  a  dielectric  plate 
streamlined  supersonic  airflow  was  designed  and  manufactured.  The  high-voltage  pulses  at  U=l- 
25  kV  with  the  help  of  a  high-voltage  cable  through  the  specially  designed  a  vacuum  cut-off 
point  is  lead  to  electrodes  fixed  into  a  dielectric  plate.  The  special  measures  for  avoidance  of  an 
electrical  breakdown  on  the  opposite  surface  of  a  dielectric  plate  outside  of  supersonic  flow  are 
undertaken.  The  plate  can  be  installed  under  different  angles  to  the  direction  of  supersonic 
airflow. 


3.4.  Automated  circuit  of  probe  measuring 

The  automated  circuit  of  probe  measuring  with  the  optical  isolator  is  designed  and 
manufactured.  The  block-scheme  of  this  circuit  is  submitted  on  Fig.  3.2. 

The  probe  current  signal  from  the  resistor  R  feeds  to  an  inlet  of  the  agreeing  amplifier  for 
amplification  up  to  value  sufficient  for  effective  coding  by  an  analog-digital  converter  (ADC).  In 
the  experiment  12-bit  4096  levels  the  analog-digital  converter  are  used.  The  time  resolution  of 
ADC  is  2  (is.  The  signals  from  the  analog-digital  converter  go  to  the  optopair  galvanic  isolation 
scheme  and  further  are  transmitted  in  a  computer  memory.  The  transmission  of  the  signals  from 


an  analog-digital  converter  into  the  computer  through  the  optical  isolator  takes  place  in  parallel 
way.  Accordingly,  in  the  circuit  12  optopairs  are  utillized  that  allows  to  not  reduce  the  operating 
speed  of  the  scheme.  The  delay  time  in  the  optopair  galvanic  isolation  scheme  does  not  exceed 
9  (is  and  the  permissible  maximum  voltage  on  its  isolation  is  7  kV.  The  potential  difference 
between  a  probe  and  antiprobe  (or  between  two  probes)  is  generated  by  the  computer.  This 
potential  difference  is  transmitted  in  a  digital  to  analog  converter  (DAC)  through  the  optopair 
galvanic  isolation  scheme  and  further  through  a  compounding  output  equipment  feeds  to  the 
probes.  The  time  resolution  of  the  DAC  is  0,5  (is,  therefore  the  basic  time  is  determined  by  the 
time  resolution  of  the  optopair  galvanic  isolation  scheme.  Thus,  the  cycle  time,  i.e.  the  time 
interval  at  which  the  circuit  is  ready  to  new  measuring  is  20  (is  at  the  intrinsic  time  of  measuring 
of  a  signal  of  a  current  t=2  (is. 

In  the  standard  version  the  scan  of  voltages  is  linear  with  sampling  of  a  step.  The  range  of 
a  voltage  variation  is  situated  in  limits  ±27  V,  the  minimum  value  of  a  step  equals  15  mV.  The 
scan  of  a  voltage  between  probes  in  principle  also  can  be  set  arbitrary,  according  to  the  program 
assigned  to  the  computer. 

The  synchronizing  circuit  allows  to  synchronize  the  pulse  of  formation  of  the  discharge 
with  a  cycle  of  measuring.  It  enables  synchronously  to  store  the  information  on  a  pulsed- 
periodical  discharge. 


Fig.  3.2.  Draft  design  of  the  set-up  for  measurement  of  voltage-current  characteristic 

of  the  double  probes. 


CHAPTER  IV 
DIAGNOSTICS  METHODS 


4.1.  Probe  application  for  diagnostics  of  the  surface  microwave  discharge 

in  supersonic  flow 

The  probe  method  permitting  to  conduct  the  measurements  of  local  plasma  parameters,  is 
represented  by  the  convenient  tool  of  examination  of  the  surface  microwave  discharges.  It  is 
concerned  with  the  fact  that  the  inhomogeneity  of  plasma  parameters  in  a  direction  both 
perpendicular,  and  parallel  surface  of  a  dielectric  surface  is  reference  feature  of  such  discharges. 

However  the  probe  method  for  diagnostics  of  such  type  of  the  discharge  was  not  applied 
yet.  Therefore  for  its  correct  application  it  is  necessary  to  analyze  problems  not  only 
interpretations  of  obtained  results,  but  also  technique  of  probe  measurements,  i.e.  select  of  an 
electric  scheme  of  a  probe  circuit,  probe  construction  and  sizes. 

At  measurements  in  a  microwave  plasma  it  is  necessary  to  take  into  account  that  for  a 
contact  method,  what  is  a  probe  method,  the  registration  of  action  of  a  very  high  frequency  field 
is  principal  moment.  Especially  it  is  important  in  a  case  of  large  pressure,  as  the  very  high 
frequency  power  which  is  necessary  for  creation  of  the  surface  microwave  discharge  grows  with 
pressure. 

Thus  two  major  factors  can  be  essential. 

At  first,  it  is  the  action  of  a  very  high  frequency  electric  fields  on  a  probe  as  a  metal 
conductor.  At  application  of  probes,  connected  with  recorded  instrumentation  by  the  radio- 
frequency  cables,  the  difficulties  arise.  They  are  bound  up  with  the  fact  that  the  conductors, 
placed  at  strong  alternating  fields,  are  the  passive  electrical  vibrators. 

Basically  it  can  result  in  to  occurrence  around  them  the  secondary  discharges  and 
essential  change  of  the  distribution  of  parameters  of  explored  plasma.  This  influence  is  most 
essential  in  a  high  frequency  region.  However  in  conditions  of  a  very  high  microwave  powers, 
taking  place  in  experiment,  the  elimination  of  this  factor  is  advisable.  Besides  in  conditions  of 
very  high  microwave  powers  it  is  desirable  to  remove  the  recorded  instrumentation  on  large 
enough  distance. 

The  second  factor  is  bound  to  action  of  a  very  high  frequency  electric  field  on  a  probe 


layer. 


Due  to  effect  of  detection  on  a  non-linearity  of  a  probe  layer  of  the  variable  electromotive 
force,  existing  between  plasma  and  ground,  a  noticeable  additional  current,  distorting  a  true 
current-voltage  characteristic  of  a  probe,  appears  in  a  probe  circuit. 

It  is  impossible  to  eliminate  this  effect  completely.  For  its  minimization  a  voltage  divider 
on  an  alternating  current  on  microwave  frequency  is  organized.  At  this  sequentially  with  a  probe 
the  filter,  which  impedance  is  much  more  than  an  impedance  of  a  probe,  is  included.  As  a  rule, 
the  quarter  wave  short-circuited  coaxial  line,  as  such  device  in  a  microwave,  is  utilized  [1]. 
However  the  probe  should  be  prolongation  of  this  line,  otherwise  the  filter  will  be  shunted. 
Therefore,  the  line  should  be  placed  in  supersonic  stream.  However  sizes  of  a  line  multiply 
exceed  the  size  of  a  probe  and  on  an  order  of  magnitude  are  comparable  with  a  dimensions  of  a 
plate.  Apparently,  it  will  lead  in  strongest  deformation  of  stream,  therefore  it  is  impossible  to  use 
such  approach. 

Let's  take  into  account,  that  the  distorted  influence  of  a  microwave  field  of  the  discharge 
takes  place  only  in  the  area  of  non-linearity  of  current-voltage  characteristic  of  a  probe.  As  a 
rule,  the  non-linearity  is  bound  up  with  an  electronic  branch  of  current-voltage  characteristic  of 
robe.  In  plasma  of  low  pressure  this  branch  is  most  informative,  allowing  in  a  number  of  cases  to 
measure  not  only  temperature  of  electrons,  but  also  electron  energy  distribution  function.  In  a 
mode  of  a  continuous  medium  the  electronic  branch  of  current-voltage  characteristic  of  probe 
can  give  the  information  only  about  conductance  of  plasma,  i.e.  information,  to  within  a  portion 
of  negative  ions,  does  not  distinguish  from  an  ionic  branch.  The  ionic  part  of  current-voltage 
characteristic  of  a  probe  in  supersonic  flow  at  large  pressures,  as  a  rule,  is  close  to  linear  [2]. 
Therefore  the  utilization  of  this  part  of  the  characteristic  allows  substantially  to  remove  a 
problem  of  distortion  of  current-voltage  characteristic. 

It  is  known  [2],  that  the  diameter  of  a  probe  being  flown  about  of  supersonic  stream 
should  be  selected  minimum,  if  only  the  mode  of  a  continuous  medium  was  executed.  However 
the  high  concentration  of  charged  particles  is  feature  of  the  surface  microwave  discharge. 
Therefore  at  a  select  of  an  ionic  branch  of  current-voltage  characteristic  the  value  of  a  useful 
signal  will  be  sufficient. 

As  the  probe  circuit  should  be  closed  one,  in  addition  to  probe  the  supporting  probe 
should  be  placed  into  explored  plasma.  Arrangement  of  the  supporting  probe  on  a  plate  will  give 
in  deformation  of  the  surface  microwave  discharge,  therefore  the  supporting  probe  needs  to  be 
disposed  in  plasma.  Taking  into  account,  that  we  were  restricted  to  an  ionic  branch,  for 
minimization  of  gasdynamic  perturbations  it  is  more  convenient  to  use  a  double  probe.  It  is 
possible  considerably  to  reduce  the  perturbations  of  investigated  system  and  recording 


instrumentation,  about  what  was  spoken  above.  For  it  an  optical  control  by  probe  power  supply 
and  an  optical  transfer  to  the  recording  instrumentation  the  information  about  a  probe  current  and 
voltage  is  used. 

Then  the  scheme  of  the  measurement  of  the  plasma  parameters  under  conditions  of  the 
surface  microwave  discharge  practically  will  not  differ  from  the  scheme  of  a  double  probe  with 
the  optical  isolator,  designed  within  the  framework  of  the  ISTC  projects  #1867p  [2]. 

The  difference  is  bound  up  with  the  fact  that  the  typical  size  of  the  surface  microwave 
discharge  in  a  direction,  perpendicular  surface  of  the  antenna  is  equal  some  millimetres.  If  you 
want  to  use  a  spherical  probe,  its  size  should  be  the  tenth  of  millimetre.  However  the  size  of  a 
leg-holder  then  should  be  the  100-th  shares  of  millimetre.  The  application  of  such  construction 
under  condition  of  supersonic  flow  practically  is  impossible.  As  a  matter  of  fact,  the  unique 
solution  in  an  explored  situation  is  the  application  of  a  double  cylindrical  probe  oriented  along  a 
surface  of  a  plate,  on  which  the  surface  microwave  discharge  is  created.  However  in  order  to  use 
the  theory  of  an  ionic  current  on  a  single  probe  under  conditions  a  double  probe,  it  is  necessary, 
that  the  interference  of  probes  was  minimum.  In  practice  it  is  reached  by  removal  of  probes  on 
distance  about  10R  (where  R  is  probe  radius),  i.e.  several  millimetres.  This  value  is  not  far  from 
the  size  of  the  discharge  in  a  transverse  direction,  therefore  such  construction  of  a  probe  is 
unacceptable.  Therefore,  it  is  possible  the  application  only  of  double  cylindrical  symmetric 
probes  oriented  along  supersonic  stream.  At  this  the  plane  of  probes  arrangement  is  parallel  the 
plane  of  a  discharge  surface. 

The  conditions  of  operation  of  a  probe  outside  of  a  boundary  layer  of  a  discharge  surface 
being  flown  about  are  close  to  conditions,  which  were  considered  in  the  ISTC  project  #1867p 
[2].  However  it  is  necessary  to  expect,  that  the  properties  of  plasma  will  change  noticeably  inside 
a  boundary  layer  at  approaching  to  a  surface.  Apart  from  change  of  a  degree  of  ionization  of 
plasma  and  temperature  its  components,  it  is  possible  change  of  a  composition  of  plasma  because 
of  an  ablation  of  a  surface.  Besides  arrangement  of  probes  along  a  stream  does  not  allow  directly 
to  utilize  results,  obtained  in  [2]  for  probes  oriented  across  of  a  stream.  The  correct  interpretation 
of  obtained  probe  characteristics  requires  the  special  additional  investigations. 

4.2.  Spectral  methods  of  measurement  of  gas  temperature 

Information  about  gas  translation  temperature  Tg  is  important  for  study  of  plasma 
processes.  The  simplest  possible  case  for  molecular  plasma  corresponds  to  local  thermodynamic 
equilibrium,  then  Tg  =  Tv  =Tr,  here  Ty  and  Tr  are  vibration  and  rotation  temperatures, 


respectively.  Ty  and  Tr  are  parameters  of  fitting  of  population  of  vibration  and  rotation  energy 
states  of  molecules  to  the  Boltzmann  function,  they  characterise  distribution  function  of 
molecules.  Both  Ty  and  Tr  can  be  extracted  from  plasma  radiation  spectra. 

In  non-equilibrium  plasmas  Tg  can  be  also  determined  on  base  of  measurement  of  Tr. 

Spectra  of  plasma  radiation  provide  information  about  rotation  temperatures  Tr  of  molecules  in 

excited  electron  states.  The  dependence  of  TR  on  intensity  of  a  rotation  structure  spectral  line  Ik, 
that  corresponds  to  a  transition  between  1-th  and  k-th  states  is  given  by  a  formula 
kR  L.  1 

Bin_jk_  =  C5a  j(j  +  i),  (4.1) 

B  lVik  Tr 

here  ke  is  the  Boltzmann  constant,  j  is  a  total  molecular  moment  of  momentum,  i  is  a  quantum 
mechanical  coefficient  of  intensity,  Vik  is  the  frequency  of  radiation  quantum,  C  is  a  constant. 

One  can  extract  TR  by  plotting  a  function  (kfi/B  )ln(I]k/iv4ik)  over  j(j+l),  and  measuring  its  slope 


tgtx: 


lR 


tga 


(4.2) 


For  a  stable  molecule  a  possibility  of  change  its  rotation  energy  at  excitation  is  negligible, 


and  rotation  distribution  functions  for  excited  and  basic  molecular  states  are  identical,  Tr=  Tr- 
Effective  exchange  of  translation  and  rotation  energies  at  collisions  results  in  equality  of  Tg  and 
Tr  of  molecules  in  basic  electron  state.  But  at  realistic  conditions  there  are  some  peculiarities 
that  must  be  paid  account  to.  If  molecules  take  part  in  plasma  chemical  reactions,  a  part  of  their 
activation  energy  can  transfer  to  rotation  energy,  that  results  in  distortion  of  initial  rotation 


distribution  function,  and  Tr^Tr.  It  is  especially  characteristic  to  chemically  active  radicals. 

Another  source  of  errors  is  bound  with  the  fact  that  spectra  yield  not  Tr  but  B  /Tr,  here  B  is  the 
rotation  constant  for  upper  state  at  radiation  transition.  At  electron  impact  excitation  of  this  state, 
rotation  energy  distribution  can  be  conserved,  but  it  corresponds  to  a  temperature 


* 


here  B°  is  the  rotation  constant  for  the  basic  state. 

Thus,  stable  molecules  with  B  ~B°  should  be  chosen  for  measurement  of  Tg.  The  second 
positive  band  of  nitrogen  is  widely  used.  Here  one  can  take  v4jk  to  be  a  constant  (the 
corresponding  error  =0.5%  is  negligible),  B  =1,826;  the  formula  for  Tr  looks  like 


(4.4) 


0,89  lg  =  C  -  — j(  j  + 1) . 

i  Tr 

In  number  of  work  it  is  suggested  to  determine  the  gas  temperature  of  nitrogen  at  high 
pressure  on  base  of  unresolved  cyan  bands  with  threshold  wavelengths  3883A0  and  3871  A0.  At 
the  present  experiments  cyan  is  generated  at  model  erosion.  Three  ways  of  temperature 
measurement  are  available.  These  methods  are  based  on:  1)  -  ratio  of  integrals  of  radiation  of  the 
two  bands;  2)  -  ratio  of  threshold  intensities  of  these  bands;  3)  -  distribution  of  rotation 
transition  intensities  of  (0;0)  band  with  threshold  wavelength  3883A0. 

4.3.  Spectral  methods  of  measurement  of  vibration  temperature 


Vibration  excitation  in  a  molecular  plasma  as  a  rule  affects  ionization  and  gas  heating. 
There  is  a  method  of  determination  of  vibration  levels  population  of  nitrogen  basic  electron  state 
X 1  Eg+.  It  is  based  on  a  measurement  of  radiation  intensities  in  bands  of  the  second  positive 
nitrogen  system.  Conditions  are  pointed  out,  at  which  the  electron  state  C  7tu  is  populated  due  to 
electron  impact  excitation  of  the  basic  X1  E„+  state.  De-excitation  of  vibration  levels  of  the  state 
C  7tu  is  bound  with  radiation,  because  its  characteristic  time  (10"  s)  is  much  less  than  that  of 
vibration  and  rotation  relaxation  and  of  diffusion. 

Let  the  vibration  levels  vx  of  the  basic  nitrogen  electron  state  are  distributed  according 
to  a  function  f  (vx ) ,  i.e. 


n(vx)  =  nof(vx)-  (4-5) 

Then  in  case  the  conditions  mentioned  above  are  met,  a  stationary  population  of  Cn  u 

vibration  levels  vc  is  given  by 

?(vx).  (4.6) 


n 


neno 


A 


GV 


'cvx 


'CVB  VX 


here  rr,  is  a  number  of  molecules  in  X1Z„+  state  with  v  Y  =  0  ;  (ov)  ■  is  a  mean  cross  section 

§  X  ’  \  /VcVx 

of  excitation  of  vc  C  7tu  levels,  ne  is  electron  concentration,  Ay  v  is  a  probability  of  radiation 

h  -j 

of  the  second  positive  nitrogen  system,  vx  are  vibration  quantum  numbers  of  B  7Cg  state. 

As  the  process  of  excitation  of  the  C?Tlu  vibration  levels  due  to  electron  impact  of  the 
basic  X'Eg+  state  is  fast  in  comparison  with  vibration  period,  these  probabilities  can  be 


considered  to  be  in  a  direct  proportion  to  the  Franc-Condon  coefficients  for  the  Tanaka  system 


C°7tu  — >  X'Lg+.  Then  the  previous  equation  takes  form 


nv  =cEq 

vc 


VCVX 


(vx). 


(4.7) 


here  C  is  a  constant  independent  on  vx .  At  computation  of  the  sum  one  can  take  only  first  5 
terms,  because  the  population  of  the  vibration  levels  is  a  steep  falling  function  of  v.  Provided 
relative  populations  of  the  first  5  levels  of  the  CM  L1  state  are  known,  one  can  calculate  relative 


populations  vx  of  the  basic  electron  states  as  a  solution  of  a  set  of  equations 

4 


n  • 

vc 


a 


q  ■  n  •• 
vcvx  vx 


(4.8) 


here  a  is  a  normalizing  coefficient.  The  relative  populations  of  the  C37tu  state  levels  can  be 
determined  experimentally  with  application  of  a  usual  method  (on  base  of  relative  intensities  of 
bands  of  the  second  positive  nitrogen  system  C37tu  — >  B37tg). 

Molecular  spectroscopy  gives  the  following  formula  for  a  radiation  intensity  of  a 
electron-vibration  band: 


I  .  =  C  Sen  ■  v4.  ■■  q  ■  -  ,  (4.9) 

VCVB  e  VC  vcvB^vcvB 

here  Se  is  a  strength  of  electron  transition,  q  are  the  Franc-Condon  factors,  v  is  a 

VCVB  VCVB 

frequency  of  transition,  C'  is  a  constant  independent  on  the  quantum  numbers  vc  and  vB  .  After 
taking  logarithm,  one  can  get 


ln(  n  •  )  =  ln(  C  )  +  In 

vc 


I  .  . 

VCVB 
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'CVB  v  C  v  B 


(4.10) 


Thus,  one  can  get  the  relative  populations  vx  and  the  vibration  temperature  of  the 
nitrogen  basic  electron  state  on  base  of  measurement  of  relative  intensities  of  the  electron- 
vibration  bands  vc. 


direction  of  supersonic  airflow  at  air  pressure  ?=40  torr,  voltage  on  the  discharge  gap  U=15  kV, 
discharge  current  i=10  A,  and  different  values  of  gas  velocity  in  supersonic  airflow.  The 
common  view  of  the  discharge  in  supersonic  airflow  on  a  surface  of  a  dielectric  body  are  given 
on  Fig. 5. 22.  It  does  not  differ  from  the  common  view  of  the  volumetric  discharge  in  supersonic 
airflow.  Without  supersonic  airflow  the  discharge  represents  the  plasma  channel  existing  on  a 
surface  of  a  dielectric  body  between  two  electrodes,  fixed  in  a  dielectric  flush  with  a  surface.  The 
supersonic  flow  leads  to  blowing  of  plasma  jets  from  each  of  electrodes  which  are  overlapped 
among  themselves  downstream.  The  length  of  each  jets  depends  on  pulse  duration,  gas  pressure, 
velocity  of  supersonic  flow,  and  energy  supplied  into  the  discharge. 

The  electron  density,  measured  by  probe  method  under  condition  transversal  pulsed- 
periodical  surface  discharge,  changed  from  ~10  cm'  at  discharge  current  i=0.5  A  up  to 
~1013  cm3  at  i=20  A. 

The  time  evolutions  of  the  gas  temperature  of  surface  pulsed-periodic  discharge  on  the 
surface  of  a  flat  plate  at  different  values  of  discharge  current  is  represented  on  Fig.  5.23.  One  can 
see,  that  at  pulse  duration  x=  1 00  ps  the  gas  temperature  increases  from  -500  K  at  discharge 
current  i=0.7  A  up  to  -3000  K  at  i=18  A. 

The  dependence  of  the  gas  heating  rate  on  discharge  current  under  condition  of  surface 
pulsed-periodic  discharge  on  the  flat  plate  is  given  on  Fig. 5. 24.  One  can  see,  that  under  our 
conditions  the  gas  heating  rate  growths  with  increasing  of  discharge  current.  This  dependence 
very  well  coordinates  with  the  gas  heating  rate  dependence  obtained  from  the  non- stationary 
kinetical  model  worked  out  by  us  [9]. 

The  vibrational  temperature  as  a  function  of  discharge  current  is  represented  on  Fig. 5. 25. 
One  can  see  that  the  vibrational  temperature  approximately  equals  5000  K  and  almost  does  not 
dependent  on  discharge  current  under  experimental  conditions. 


Front  view  Side  view 

Fig. 5. 22.  The  common  view  of  the  pulsed-periodic  discharge  in  supersonic  airflow  on  a  surface 
of  flat  plane  at  different  values  of  air  pressure  in  the  high-pressure  system  p,  atm:  1-1;  2-3;  3-6. 


Fig. 5. 23.  The  gas  temperature  of  surface  pulsed-periodic  discharge  on  the  flat  plate  at  different 


values  of  discharge  current  i,  A:  1  -  0.7;  2  -  1.5;  3  -  3.0;  4  -  6.0;  5  -  9.0;  6-18. 


Fig. 5. 24.  Dependence  of  the  gas  heating  rate  on  discharge  current  under  condition  of  surface 

pulsed-periodic  discharge  on  the  flat  plate. 


Fig. 5. 25.  Dependence  of  the  vibrational  temperature  on  discharge  current  under  condition 
of  surface  pulsed-periodic  discharge  on  the  flat  plate. 


CONCLUSIONS 


1.  The  literature  review  on  the  subject  of  the  Project  is  carried  out.  The  different  regimes  of 
action  on  laminar  -  turbulent  transition  and  turbulent  boundary  layer  resulting  in  to  reduction 
of  friction  at  subsonic  and  supersonic  motion  of  flight  vehicles  in  the  Earth's  atmosphere  are 
analysed.  The  different  ways  of  additional  energy  supply  into  a  turbulent  boundary  layer  are 
examined.  The  different  modes  of  the  gas  discharges  creation  into  a  boundary  layer  are 
considered.  The  methods  of  plasma  creation  on  the  bodies  surfaces  with  the  help  of  the 
surface  microwave  discharges  are  examined. 

2.  The  mathematical  model  of  the  discharge  grounded  on  the  joint  solution  of  the  balance 
equations  of  charged  particles  and  Maxwell  equations  for  an  electromagnetic  field  are 
worked  out.  The  theoretical  analysis  ensures  following  results:  the  calculation  of  the 
dispersion  curves  of  a  surface  wave  taking  into  account  inhomogeneity  of  a  plasma  layer 
(self-consistent  solution  with  the  ionised  non-linearity);  the  calculation  of  a  field  penetration 
depth  in  plasma  and  thickness  of  a  plasma  layer  created  by  both  direct  and  reflected  surface 
waves,  in  dependence  on  wave  amplitude  and  frequency,  and  also  gas  pressure  and  airflow 
velocity;  the  dependence  of  longitudinal  length  of  the  discharge  supported  by  a  surface  wave 
from  power,  transferred  by  a  microwave;  the  analysis  of  requirements  of  electromagnetic 
waves  reflection  from  antenna  butt-end  at  the  different  shape  of  the  latter;  the  analysis  of  two 
dimensional  field  structure  in  a  plane  of  the  antenna. 

3.  The  modernization  of  the  discharge  chamber  available  in  laboratory  is  carried  out  with  the 
purpose  of  creation  of  the  surface  microwave  and  pulsed-periodic  discharges  in  supersonic 
airflow  at  gas  pressure  p=l-760  torr.  The  system  of  formation  of  supersonic  gas  flow  and  the 
special  system  of  input  of  a  microwave  energy  into  discharge  chamber  and  transformations 
of  the  microwave  energy  into  the  energy  of  a  surface  wave  is  installed  on  the  discharge 
chamber.  The  method  of  probe  measurements  of  the  charged  particles  concentration  under 
conditions  of  surface  microwave  discharge  plasma  was  worked  out.  It  is  shown  that  the 
perturbations  of  researched  system  and  distortion  of  the  probe  characteristics  can 
considerably  be  reduced  by  excluding  the  radio-frequency  cables  from  constructions  of 
devise  and  by  supplied  optical  management  of  a  probe  feeding  and  transfer  of  the 
information  about  a  probe  current  and  voltage  into  the  recording  equipment.  The  automated 
circuit  of  probe  measuring  with  the  optical  isolator  is  designed  and  manufactured.  In  the 
experiment  12-bit  4096  levels  the  analog-digital  converter  are  used. 


4.  The  surface  microwave  discharge  on  a  dielectric  antenna  with  length  L=15  cm  of  rectangular 
section  S=l?2  cm2  is  created  at  air  pressure  ?=10~3-103  torr.  The  magnetron  generator  had  the 
following  characteristics:  the  wavelength  A,=2,4  cm;  the  pulsed  microwave  power  W=10- 
300  kW;  the  repetition  frequency  f=l-100Hz;  the  pulse  duration  x=  1  -f  300  ps;  the  pulse 
period-to-pulse  duration  ratio  Q=1000.  At  air  pressure  ?>1  torr  the  discharge  represents  thin 
(thickness  Ah=0, 1-1,0  mm)  plasma  layer  uniformly  coating  all  antenna  surface.  The 
examination  of  dynamics  of  the  discharge  has  shown  that  velocity  of  propagation  of  the 
discharge  and  its  longitudinal  length  depend  on  microwave  pulse  duration  and  input  power. 

n 

Thus,  on  an  initial  stage  of  the  discharge  the  velocity  of  its  distribution  reaches  v~10  cm/s, 
that  are  much  greater  than  the  airflow  velocity  (vf=5-104  cm/s)  under  experimental 
conditions.  The  discharge  propagation  velocity  decreases  up  to  v~104  cm/s  to  the  end  of  the 
microwave  pulse  with  duration  x=100  (is.  The  threshold  characteristic  of  the  surface 
microwave  discharge  appearance  are  determined.  It  is  shown,  that  at  air  pressure  ?=0,1- 
10  torr  the  power  necessary  for  creation  of  the  surface  microwave  discharge  is  minimum  and 
equals  ~10  kW. 

5.  In  plasma  of  the  surface  microwave  discharge  the  electron  density,  gas  and  vibrational 
temperature  are  measured  at  different  values  of  gas  pressure  and  microwave  power.  The 
electron  density  equals  610  cm"  .  At  growth  of  power  the  gas  temperature  is  increased 
from  -500  K  at  W=35  kW  up  to  -1700  K  at  W=175  kW,  whereas  vibrational  temperature 
remains  practically  constant  under  these  conditions,  insignificantly  decreasing  with  increase 
of  microwave  power.  At  this  the  maximal  gas  heating  is  observed  in  a  place  of  excitation  of  a 
surface  microwave  discharge  and  gas  temperature  decreases  by  the  end  of  the  discharge.  On 
the  initial  stage  of  existence  of  the  surface  microwave  discharge  the  fast  gas  heating  with  rate 
dTg/dt=50  K/|is  is  observed.  It  was  shown,  that  the  mechanism  connected  to  effective 
excitation  of  electron-exited  states  of  nitrogen  molecules  at  large  values  of  the  reduced 
electrical  field  E/ri>100  Td=10~15  Vcm2  and  their  subsequent  quenching  is  responsible  for 
the  fast  heating.  At  this  the  part  of  excitation  energy  of  these  states  is  transferred  in  the  heat 
of  air.  As  our  estimations  show,  only  the  quenching  of  electron-exited  long-living  states  of 
nitrogen  molecules,  which  are  effectively  created  in  conditions  of  a  surface  microwave 
discharge  in  air,  provides  the  gas  heating  rate. 

6 .  The  surface  microwave  discharges  in  supersonic  airflow  with  Mach  number  of  flow  ?  =2  are 
created.  It  was  shown,  that  the  supersonic  flow  of  air  at  Mach  number  ?  =2  does  not 
influence  on  the  common  view  of  the  surface  discharge  on  the  dielectric  antenna  of 


rectangular  section  and  on  the  microwave  power  which  is  necessary  for  creation  of  the 
surface  microwave  discharge.  Thus  the  threshold  of  appearance  of  a  surface  microwave 
discharge  does  not  depend  on  whether  at  first  the  discharge  on  a  surface  of  a  dielectric  body 
is  created  and  then  the  supersonic  airflow  is  included  or,  on  the  contrary,  at  first  the 
supersonic  airflow  is  included,  and  only  then  the  surface  microwave  discharge  is  created. 

7.  The  system  of  creation  of  a  pulsed-periodic  discharge  on  a  surface  of  a  dielectric  plate 
streamlined  supersonic  airflow  was  designed  and  manufactured.  The  high-voltage  pulses  at 
U=l-25  kV  with  the  help  of  a  high-voltage  cable  through  the  specially  designed  a  vacuum 
cut-off  point  is  lead  to  electrodes  fixed  into  a  dielectric  plate.  The  special  measures  for 
avoidance  of  an  electrical  breakdown  on  the  opposite  surface  of  a  dielectric  plate  outside  of 
supersonic  flow  are  undertaken.  The  plate  can  be  installed  under  different  angles  to  the 
direction  of  supersonic  airflow. 

8 .  The  surface  pulsed-periodic  discharge  in  supersonic  airflow  was  created  on  a  dielectric  plate 
placed  at  the  angle  of  approximately  twenty  degrees  to  the  direction  of  supersonic  airflow  at 
air  pressure  ?=40  torr,  voltage  on  the  discharge  gap  U=15  kV,  and  discharge  current  i=10  A. 
The  common  view  of  the  discharge  in  supersonic  airflow  on  a  surface  of  a  dielectric  body 
does  not  differ  from  the  common  view  of  the  volumetric  discharge  in  supersonic  airflow. 
Without  supersonic  airflow  the  discharge  represents  the  plasma  channel  existing  on  a  surface 
of  a  dielectric  body  between  two  electrodes,  fixed  in  a  dielectric  flush  with  a  surface.  The 
supersonic  flow  leads  to  blowing  of  plasma  jets  from  each  of  electrodes  which  are 
overlapped  among  themselves  downstream.  The  length  of  each  jets  depends  on  pulse 
duration,  gas  pressure,  velocity  of  supersonic  flow,  and  energy  deposition  in  the  discharge. 
The  electron  density,  measured  by  probe  method  under  condition  transversal  pulsed- 
periodical  surface  discharge,  changed  from  ~10  cm"  at  discharge  current  i=0.5  A  up  to 
~  10  cm"  at  i=20  A.  The  time  evolutions  of  the  gas  temperature  of  surface  pulsed-periodic 
discharge  on  the  surface  of  a  flat  plate  at  different  values  of  discharge  current  is  obtained.  It 
is  shown,  that  under  experimental  conditions  the  gas  heating  rate  growths  with  increasing  of 
discharge  current.  The  vibrational  temperature  approximately  equals  5000  K. 
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direction  of  supersonic  airflow  at  air  pressure  ?=40  torr,  voltage  on  the  discharge  gap  U=15  kV, 
discharge  current  i=10  A,  and  different  values  of  gas  velocity  in  supersonic  airflow.  The 
common  view  of  the  discharge  in  supersonic  airflow  on  a  surface  of  a  dielectric  body  are  given 
on  Fig. 5. 22.  It  does  not  differ  from  the  common  view  of  the  volumetric  discharge  in  supersonic 
airflow.  Without  supersonic  airflow  the  discharge  represents  the  plasma  channel  existing  on  a 
surface  of  a  dielectric  body  between  two  electrodes,  fixed  in  a  dielectric  flush  with  a  surface.  The 
supersonic  flow  leads  to  blowing  of  plasma  jets  from  each  of  electrodes  which  are  overlapped 
among  themselves  downstream.  The  length  of  each  jets  depends  on  pulse  duration,  gas  pressure, 
velocity  of  supersonic  flow,  and  energy  supplied  into  the  discharge. 

The  electron  density,  measured  by  probe  method  under  condition  transversal  pulsed- 
periodical  surface  discharge,  changed  from  ~10  cm'  at  discharge  current  i=0.5  A  up  to 
~1013  cm3  at  i=20  A. 

The  time  evolutions  of  the  gas  temperature  of  surface  pulsed-periodic  discharge  on  the 
surface  of  a  flat  plate  at  different  values  of  discharge  current  is  represented  on  Fig.  5.23.  One  can 
see,  that  at  pulse  duration  x=  1 00  ps  the  gas  temperature  increases  from  -500  K  at  discharge 
current  i=0.7  A  up  to  -3000  K  at  i=18  A. 

The  dependence  of  the  gas  heating  rate  on  discharge  current  under  condition  of  surface 
pulsed-periodic  discharge  on  the  flat  plate  is  given  on  Fig. 5. 24.  One  can  see,  that  under  our 
conditions  the  gas  heating  rate  growths  with  increasing  of  discharge  current.  This  dependence 
very  well  coordinates  with  the  gas  heating  rate  dependence  obtained  from  the  non- stationary 
kinetical  model  worked  out  by  us  [9]. 

The  vibrational  temperature  as  a  function  of  discharge  current  is  represented  on  Fig. 5. 25. 
One  can  see  that  the  vibrational  temperature  approximately  equals  5000  K  and  almost  does  not 
dependent  on  discharge  current  under  experimental  conditions. 


Front  view  Side  view 

Fig. 5. 22.  The  common  view  of  the  pulsed-periodic  discharge  in  supersonic  airflow  on  a  surface 
of  flat  plane  at  different  values  of  air  pressure  in  the  high-pressure  system  p,  atm:  1-1;  2-3;  3-6. 


Fig. 5. 23.  The  gas  temperature  of  surface  pulsed-periodic  discharge  on  the  flat  plate  at  different 


values  of  discharge  current  i,  A:  1  -  0.7;  2  -  1.5;  3  -  3.0;  4  -  6.0;  5  -  9.0;  6-18. 


Fig. 5. 24.  Dependence  of  the  gas  heating  rate  on  discharge  current  under  condition  of  surface 

pulsed-periodic  discharge  on  the  flat  plate. 


Fig. 5. 25.  Dependence  of  the  vibrational  temperature  on  discharge  current  under  condition 
of  surface  pulsed-periodic  discharge  on  the  flat  plate. 


CONCLUSIONS 


1.  The  literature  review  on  the  subject  of  the  Project  is  carried  out.  The  different  regimes  of 
action  on  laminar  -  turbulent  transition  and  turbulent  boundary  layer  resulting  in  to  reduction 
of  friction  at  subsonic  and  supersonic  motion  of  flight  vehicles  in  the  Earth's  atmosphere  are 
analysed.  The  different  ways  of  additional  energy  supply  into  a  turbulent  boundary  layer  are 
examined.  The  different  modes  of  the  gas  discharges  creation  into  a  boundary  layer  are 
considered.  The  methods  of  plasma  creation  on  the  bodies  surfaces  with  the  help  of  the 
surface  microwave  discharges  are  examined. 

2.  The  mathematical  model  of  the  discharge  grounded  on  the  joint  solution  of  the  balance 
equations  of  charged  particles  and  Maxwell  equations  for  an  electromagnetic  field  are 
worked  out.  The  theoretical  analysis  ensures  following  results:  the  calculation  of  the 
dispersion  curves  of  a  surface  wave  taking  into  account  inhomogeneity  of  a  plasma  layer 
(self-consistent  solution  with  the  ionised  non-linearity);  the  calculation  of  a  field  penetration 
depth  in  plasma  and  thickness  of  a  plasma  layer  created  by  both  direct  and  reflected  surface 
waves,  in  dependence  on  wave  amplitude  and  frequency,  and  also  gas  pressure  and  airflow 
velocity;  the  dependence  of  longitudinal  length  of  the  discharge  supported  by  a  surface  wave 
from  power,  transferred  by  a  microwave;  the  analysis  of  requirements  of  electromagnetic 
waves  reflection  from  antenna  butt-end  at  the  different  shape  of  the  latter;  the  analysis  of  two 
dimensional  field  structure  in  a  plane  of  the  antenna. 

3.  The  modernization  of  the  discharge  chamber  available  in  laboratory  is  carried  out  with  the 
purpose  of  creation  of  the  surface  microwave  and  pulsed-periodic  discharges  in  supersonic 
airflow  at  gas  pressure  p=l-760  torr.  The  system  of  formation  of  supersonic  gas  flow  and  the 
special  system  of  input  of  a  microwave  energy  into  discharge  chamber  and  transformations 
of  the  microwave  energy  into  the  energy  of  a  surface  wave  is  installed  on  the  discharge 
chamber.  The  method  of  probe  measurements  of  the  charged  particles  concentration  under 
conditions  of  surface  microwave  discharge  plasma  was  worked  out.  It  is  shown  that  the 
perturbations  of  researched  system  and  distortion  of  the  probe  characteristics  can 
considerably  be  reduced  by  excluding  the  radio-frequency  cables  from  constructions  of 
devise  and  by  supplied  optical  management  of  a  probe  feeding  and  transfer  of  the 
information  about  a  probe  current  and  voltage  into  the  recording  equipment.  The  automated 
circuit  of  probe  measuring  with  the  optical  isolator  is  designed  and  manufactured.  In  the 
experiment  12-bit  4096  levels  the  analog-digital  converter  are  used. 


4.  The  surface  microwave  discharge  on  a  dielectric  antenna  with  length  L=15  cm  of  rectangular 
section  S=l?2  cm2  is  created  at  air  pressure  ?=10~3-103  torr.  The  magnetron  generator  had  the 
following  characteristics:  the  wavelength  A,=2,4  cm;  the  pulsed  microwave  power  W=10- 
300  kW;  the  repetition  frequency  f=l-100Hz;  the  pulse  duration  x=  1  -f  300  ps;  the  pulse 
period-to-pulse  duration  ratio  Q=1000.  At  air  pressure  ?>1  torr  the  discharge  represents  thin 
(thickness  Ah=0, 1-1,0  mm)  plasma  layer  uniformly  coating  all  antenna  surface.  The 
examination  of  dynamics  of  the  discharge  has  shown  that  velocity  of  propagation  of  the 
discharge  and  its  longitudinal  length  depend  on  microwave  pulse  duration  and  input  power. 

n 

Thus,  on  an  initial  stage  of  the  discharge  the  velocity  of  its  distribution  reaches  v~10  cm/s, 
that  are  much  greater  than  the  airflow  velocity  (vf=5-104  cm/s)  under  experimental 
conditions.  The  discharge  propagation  velocity  decreases  up  to  v~104  cm/s  to  the  end  of  the 
microwave  pulse  with  duration  x=100  (is.  The  threshold  characteristic  of  the  surface 
microwave  discharge  appearance  are  determined.  It  is  shown,  that  at  air  pressure  ?=0,1- 
10  torr  the  power  necessary  for  creation  of  the  surface  microwave  discharge  is  minimum  and 
equals  ~10  kW. 

5.  In  plasma  of  the  surface  microwave  discharge  the  electron  density,  gas  and  vibrational 
temperature  are  measured  at  different  values  of  gas  pressure  and  microwave  power.  The 
electron  density  equals  610  cm"  .  At  growth  of  power  the  gas  temperature  is  increased 
from  -500  K  at  W=35  kW  up  to  -1700  K  at  W=175  kW,  whereas  vibrational  temperature 
remains  practically  constant  under  these  conditions,  insignificantly  decreasing  with  increase 
of  microwave  power.  At  this  the  maximal  gas  heating  is  observed  in  a  place  of  excitation  of  a 
surface  microwave  discharge  and  gas  temperature  decreases  by  the  end  of  the  discharge.  On 
the  initial  stage  of  existence  of  the  surface  microwave  discharge  the  fast  gas  heating  with  rate 
dTg/dt=50  K/|is  is  observed.  It  was  shown,  that  the  mechanism  connected  to  effective 
excitation  of  electron-exited  states  of  nitrogen  molecules  at  large  values  of  the  reduced 
electrical  field  E/ri>100  Td=10~15  Vcm2  and  their  subsequent  quenching  is  responsible  for 
the  fast  heating.  At  this  the  part  of  excitation  energy  of  these  states  is  transferred  in  the  heat 
of  air.  As  our  estimations  show,  only  the  quenching  of  electron-exited  long-living  states  of 
nitrogen  molecules,  which  are  effectively  created  in  conditions  of  a  surface  microwave 
discharge  in  air,  provides  the  gas  heating  rate. 

6 .  The  surface  microwave  discharges  in  supersonic  airflow  with  Mach  number  of  flow  ?  =2  are 
created.  It  was  shown,  that  the  supersonic  flow  of  air  at  Mach  number  ?  =2  does  not 
influence  on  the  common  view  of  the  surface  discharge  on  the  dielectric  antenna  of 


rectangular  section  and  on  the  microwave  power  which  is  necessary  for  creation  of  the 
surface  microwave  discharge.  Thus  the  threshold  of  appearance  of  a  surface  microwave 
discharge  does  not  depend  on  whether  at  first  the  discharge  on  a  surface  of  a  dielectric  body 
is  created  and  then  the  supersonic  airflow  is  included  or,  on  the  contrary,  at  first  the 
supersonic  airflow  is  included,  and  only  then  the  surface  microwave  discharge  is  created. 

7.  The  system  of  creation  of  a  pulsed-periodic  discharge  on  a  surface  of  a  dielectric  plate 
streamlined  supersonic  airflow  was  designed  and  manufactured.  The  high-voltage  pulses  at 
U=l-25  kV  with  the  help  of  a  high-voltage  cable  through  the  specially  designed  a  vacuum 
cut-off  point  is  lead  to  electrodes  fixed  into  a  dielectric  plate.  The  special  measures  for 
avoidance  of  an  electrical  breakdown  on  the  opposite  surface  of  a  dielectric  plate  outside  of 
supersonic  flow  are  undertaken.  The  plate  can  be  installed  under  different  angles  to  the 
direction  of  supersonic  airflow. 

8 .  The  surface  pulsed-periodic  discharge  in  supersonic  airflow  was  created  on  a  dielectric  plate 
placed  at  the  angle  of  approximately  twenty  degrees  to  the  direction  of  supersonic  airflow  at 
air  pressure  ?=40  torr,  voltage  on  the  discharge  gap  U=15  kV,  and  discharge  current  i=10  A. 
The  common  view  of  the  discharge  in  supersonic  airflow  on  a  surface  of  a  dielectric  body 
does  not  differ  from  the  common  view  of  the  volumetric  discharge  in  supersonic  airflow. 
Without  supersonic  airflow  the  discharge  represents  the  plasma  channel  existing  on  a  surface 
of  a  dielectric  body  between  two  electrodes,  fixed  in  a  dielectric  flush  with  a  surface.  The 
supersonic  flow  leads  to  blowing  of  plasma  jets  from  each  of  electrodes  which  are 
overlapped  among  themselves  downstream.  The  length  of  each  jets  depends  on  pulse 
duration,  gas  pressure,  velocity  of  supersonic  flow,  and  energy  deposition  in  the  discharge. 
The  electron  density,  measured  by  probe  method  under  condition  transversal  pulsed- 
periodical  surface  discharge,  changed  from  ~10  cm"  at  discharge  current  i=0.5  A  up  to 
~  10  cm"  at  i=20  A.  The  time  evolutions  of  the  gas  temperature  of  surface  pulsed-periodic 
discharge  on  the  surface  of  a  flat  plate  at  different  values  of  discharge  current  is  obtained.  It 
is  shown,  that  under  experimental  conditions  the  gas  heating  rate  growths  with  increasing  of 
discharge  current.  The  vibrational  temperature  approximately  equals  5000  K. 
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